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Abstract. We study the out-of-equilibrium dynamics of the spherical 
ferromagnet after a quench to its critical temperature. We calculate correlation 
and response functions for spin observables which probe lengthscales much larger 
than the lattice spacing but smaller than the system size, and find that the 
asymptotic fluctuation-dissipation ratio (FDR) X°° is the same as for local 
observables. This is consistent with our earlier results for the Ising model in 
dimension d = 1 and d = 2. We also check that bond observables, both local 
and long-range, give the same asymptotic FDR. In the second part of the paper 
the analysis is extended to global observables, which probe correlations among 
all N spins. Here non-Gaussian fluctuations arising from the spherical constraint 
need to be accounted for, and we develop a systematic expansion in l/y/N to 
do this. Applying this to the global bond observable, i.e. the energy, we find 
that non-Gaussian corrections change its FDR to a nontrivial value which we 
calculate exactly for all dimensions d > 2. Finally, we consider quenches from 
magnetized initial states. Here even the FDR for the global spin observable, i.e. 
the magnetization, is nontrivial. It differs from the one for unmagnetized states 
even in d > 4, signalling the appearance of a distinct dynamical universality class 
of magnetized critical coarsening. For lower d the FDR is irrational even to first 
order in 4 — d and d — 2, the latter in contrast to recent results for the n- vector 
model. 



1. Introduction 

A key insight of statistical mechanics is that equilibrium states can be accurately 
described in terms of only a small number of thermodynamic variables, such as 
temperature and pressure. For non-equilibrium systems such as glasses no similar 
simplification exists a priori; the whole past history of a sample is in principle required 
to specify its state at a given time. This complexity makes theoretical analysis 
awkward, and one is led instead to look for a description of non-equilibrium states in 
terms of a few effective thermodynamic parameters. Much work in recent years has 
focussed on one such parameter, the effective temperature. This can be defined on 
the basis of fluctuation-dissipation (FD) relations between correlation and response 
functions and has proved to be very fruitful in mean field systems Pj |2| ■ 

The use of FD relations to quantify the out-of-equilibrium dynamics in glassy 
systems is motivated by the occurrence of aging the time scale of response 
to an external perturbation increases with the age (time since preparation) t„ of 

t Email alessia.annibale@kcl.ac.uk 
X Email peter.sollich@kcl.ac.uk 



Fluctuation-dissipation relations in the non- equilibrium spherical ferromagnet 



2 



the system. As a consequence, time translational invariance and the equihbrium 
fluctuation-dissipation theorem ^ (FDT) relating correlation and response functions 
break down. To quantify this, one considers the correlation function of two generic 
observable A and B oi a, system, defined as 

C{t, t^) = {A{t)B{t^)) - {A{t)) {B{t^)) (1.1) 

The associated (impulse) response function can be defined as 

5{A{t)) 



R{t, tw) 



Shsit^) 



and gives the linear response of A at time i to a small impulse in the field hs conjugate 
to B at the earlier time (The latter is normally thought of as a "waiting time" since 
preparation of the system at time 0.) Equivalently one can work with the susceptibility 

X(i,tw)= / dt'R{t,t') (1.2) 



which encodes the response of A{t) to a small step hB{t) — hsQit — tw) in the 
field starting at t^. In equilibrium, FDT implies that —dt„x{t,t^) — R(t,t^) ~ 
T^^dt„C{t, tw). Out of equilibrium, the violation of FDT can be measured by an FD 
ratio (FDR), X, defined through OE] 

-9*„x(t,iw) = i?(t,tw) = ^^^|^9*„C(t,tw) (1.3) 

This implies that X can be read off from the slope ~X/T of a parametric FD plot 
showing X vs C, at fixed t and with t^ as the curve parameter. This remains the case 
if both axes are normalized by the equal-time variance of A^ C{t, t), a procedure which 
is helpful in fixing the scale of the plot in situations where C{t, t) varies significantly 
with time jTUHj. In equilibrium, the FD plot is a straight line with slope —1/T. 

In mean-field spin glasses [2 El El one finds that FD plots of autocorrelation and 
response of local spins and similar observables approach a limiting shape for large t. 
This is typically composed of two straight line segments. In the first of these one 
finds X — \, corresponding to quasi-equilibrium dynamics for time differences t — ty^ 
that do not grow with the age of the system. The second line segment has X < 1 
and reflects the dynamics on aging timescales, i.e. time differences growing (in the 
simplest case linearly) with t^. One can use this to deflne a non-equilibrium effective 
temperature Teg = T/X, which has been shown to have many of the properties of a 
thermodynamic temperature j^ElEl- 

How well this physically very attractive mean-field scenario transfers to more 
realistic non-equilibrium systems with short-range interactions has been a matter of 
intense research recently A useful class of systems for studying this question in 
detail is provided by ferromagnets quenched from high temperature to the critical 
temperature (see e.g. Refs. (51 El and the recent review [H] ) or below. The system 
then coarsens - by the growth of domains with the equilibrium magnetization, for 
T < Tc - and exhibits aging; in an infinite system equilibrium is never reached. The 
aging is clearly related to the growth of a lengthscale ^2 (domain size for T < Tc, or 
correlation length for T = Tc), and this makes ferromagnets attractive "laboratory" 
systems for understanding general properties of non-equilibrium dynamics. They are 
of course not completely generic; compared to e.g. glasses they lack features such as 
thermal activation over energetic or entropic barriers. 
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We focus in this paper mostly on ferromagnets quenched to Tc, i.e. on critical 
coarsening dynamics. Some care is needed in this case with the interpretation of 
limiting FD plots: while in the mean-field situation X becomes at long times a function 
of C only, as implied by the existence of a nontrivial limit plot, in critical coarsening X 
approaches a function oit/t^ In the interesting regime where t/ty^ is finite but > 1, 
time differences t — tvj = C'(iw) are then large, and e.g. spin autocorrelation functions 
have decayed to a small part of their initial value. In the limit — oo the FD plot 
then assumes a pseudo-equilibrium shape, with all nontrivial detail compressed into a 
vanishing region around C — 0. 

The fact that the FDR is a smooth function of t/ty^ makes the interpretation of 
T/X as an effective temperature less obvious than in mean- field spin glasses, where 
T/X is constant within each time sector (t — t^, = 0(1) vs t — growing with t^)- 
To eliminate the time-dependence one can consider the limit of times that are both 
large and well-separated. This defines an asymptotic FDR 

X°° = lim lim X{t,ty,) (1.4) 

tw — *00 t — >00 

An important property of this quantity is that it should be universal ^ lll| in the 
sense that its value is the same for different systems falling into the same universality 
class of critical non-equilibrium dynamics. This makes a study of X°° interesting in 
its own right, even without an interpretation in terms of effective temperatures. 

If one nevertheless wants to pursue such an interpretation, the resulting value 
of the effective temperature T/X°° should be the same for all (or at least a large 
class of) observables A. The observable-dependence oi X°° therefore becomes a key 
question [3 |H1 1101 113| . Conventionally, much work on non-equilibrium ferromagnets 
has focussed on the local spin autocorrelation function and associated response. An 
obvious alternative is the long-wavelength analogue, i.e. the correlation function of the 
fluctuating magnetization. Exact calculations for the Ising chain |1U[ [T^ as well as 
numerical simulations jlOi, il5j in dimension d = 2 show that the resulting global X°° 
is always identical to the local version. This local-global correspondence, which can 
also be obtained by field-theoretic arguments jlllllfillTTj . arises physically because the 
long wavelength Fourier components of the spins are slowest to relax and dominate 
the long-time behaviour of both local and global quantities. 

The local-global correspondence docs of course not address the full range of 
observable-dependence of the asymptotic FDR; one might ask about other observables 
which are linear combinations not of spins but for example products of interacting 
spins. In the critical Ising model in d = 2, numerical simulations El suggest 
that even these give the same X°°, so that an interpretation of T/X°" in terms of 
an effective temperature appears plausible. One of the motivations for the current 
study was to verify whether this observable-independence of X°° across different types 
of observables holds in an exactly solvable model, the spherical ferromagnet 18, 1^. 
In addition, we will study what effect different initial conditions have on X°°. This 
is motivated by our recent study of Ising models in the classical regime of large d, 
where critical fluctuations are irrelevant |2(J| . It turned out that magnetized initial 
do states produce a different value of X°°, so that critical coarsening in the presence 
of a nonzero magnetization is in a new dynamical universality class even though the 
magnetization does decay to zero at long times. 

We begin in Sec. |21 with a brief review of the standard setup for the dynamics of 
the spherical model, as used in e.g. jH]. Fluctuations in an effective Lagrange multiplier 
enforcing the spherical constraint are neglected, leading to a theory where all spins are 
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Gaussian random variables. In Sec. 13 this is applied to various observables of finite 
range, by which we mean correlations and responses probing lengthscales that can be 
large but remain small compared to the overall system size. For spin observables we 
show that the expected equality of X'^ between local and long-range quantities holds 
fSec. lB.lll . We check observable-independence of X°° further by considering bond and 
spin product observables, in Sec . 13 . 21 and 13.31 resoectivelv. 

The major part of the paper is then devoted to a study of FDRs for global 
observables, with a focus on the energy, i.e. the global bond observable. Because 
of the weak infinite-range interaction generated by the spherical constraint, such 
observables behave differently from their long-range analogues in the spherical model. 
Calculations of correlation and response functions are technically substantially more 
difficult because Lagrange multiplier fluctuations can no longer be neglected. To 
account for them we construct in Sec. 21 8. systematic expansion of the dynamically 
evolving spins in A^~^/^. This allows us to calculate the leading non-Gaussian 
corrections that we need for global correlations, as shown for the case of the energy 
in Sec. ]E\ After a brief digression to equilibrium dynamics, we evaluate the resulting 
expressions in SecElfor d above the critical dimension dc — 4, and in Sec.[3for d < 4. 
Importantly, we will find that in the latter case the asymptotic FDR is different 
from those for finite-range observable. This means that an effective temperature 
interpretation of X°° is possible at best in a very restricted sense. However, we 
will find that our results are in agreement with recent renormalization group (RG) 
calculations near d = 4 ^3] in the 0{n (X))-model. This suggests that the non- 
Gaussian effects captured in global observables are important for linking the spherical 
model to more realistic systems with only short-range interactions. Finally, we turn in 
Sec.|Slto critical coarsening starting from magnetized initial conditions. Here already 
the global spin observable is affected by non-Gaussian corrections. Once these are 
accounted for, we find X°° = 4/5 for d > 4 as in the Ising case [201 • For d < 4 we 
provide the first exact values of the asymptotic FDR in the presence of a nonzero 
magnetization; these turn out to be highly nontrivial even to first order in 4 — d and 
d — 2. Our results are summarized and discussed in Sec. [HI Technical details are 
relegated to two appendices. 

2. Langevin dynamics and Gaussian theory 

We consider the standard spherical model Hamiltonian 

H=lY.(^^-S,f (2.1) 

iv) 

The sum runs over all nearest neighbour (n.n.) pairs on a d-dimensional (hyper-)cubic 
lattice; the lattice constant is taken as the unit length. At each of the N lattice sites 
there is a real- valued spin Si. The spherical constraint Sf ^ N is imposed, which 
can be motivated by analogy with Ising spins Si = ±1 |18|. 
The Langevin dynamics for this model can be written as 

with ^i Gaussian white noise with zero mean and covariance {S.i{t)^j{t')) = 2T6ijS{t — 
t'). The last term in (|2.2|l . i.e. the sum over k, enforces the spherical constraint at 
all times by removing the component of the velocity vector {dtSi, . . . ,dtSN) along 
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(5*1, . . . , S'at). We use here the Stratonovic convention for products Uke SkS,k- This 
allows the ordinary rules of calculus to be used when evaluating derivatives such as 
dtSf. PhysicaUy it corresponds to the intuitively reasonable scenario where the noise 
ifi is regarded as a smooth random process but with a correlation time much shorter 
than any other dynamical timescale. 

The prefactor of Si in the last term of 1)2. 2|l . being an average of N 
contributions, will have fluctuations of 0(7V~^/^). Conventionally one ignores these 
and approximates the equation of motion as 

B H 

^tS^ = -— + C^~z{t)S, (2.3) 
aSi 

where z(t) can be viewed as an effective time-dependent Lagrange multiplier 
implementing the spherical constraint. This approximation works for local quantities, 
but as we will see can give incorrect results when one considers e.g. fluctuations of 
the magnetization or the energy, which involve correlations across the entire system. 
One can see directly that 1)2. 3|l is an approximation from the fact that it corresponds 
to Langevin dynamics with the effective Hamiltonian H + ■^z{t)^- Sf. Since the 
latter is time-dependent, this dynamics does not satisfy detailed balance. It is 
simple to check, on the other hand, that the original equation of motion (|2.2|l does 
satisfy detailed balance and leads to the correct equilibrium distribution Poq({*S'i}) oc 
exp{—(3H)6{J2 Sf — N) where jS — 1/T is the inverse temperature as usual. 

The key advantage of the approximation 1)2.3(1 is, of course, that the spins are 
Gaussian random variables at all times as long as the initial condition is of this form. 
Explicitly, if we define a matrix ft with Qij = — 1 for n.n. sites i,j and Qu = 2d, the 
Gaussian equation of motion is 

dtS, = - ^ - z(t)S, + 6 (2.4) 
j 

We review briefly how this is solved (see e.g. and references therein), since these 
results form the basis for all later developments. In terms of the Fourier components 
Sq = Si exp(— jq • Ti) of the spins, equation ()2.4)) reads 

9t5q = -(c^q + z(i))S'q+eq (2-5) 

where ujq — 2X]f=i(l ~ cosga); we mostly write just to. The Fourier mode response 
function can be read off as 



i?q(i, iw) = exp (^-^u;{t - t^) - J* dt' zit')^ = J ^e-'^(*-*-) 



(2.6) 



where 

g{t) = exp ( 2 / dt' z{t') ] (2.7) 







In terms of this, the time-dependence of the Sc^ becomes 

S'q(t) = i?q(t, O)S'q(O) + / dt' Rq{t, t')Cq{t') (2., 







The equal-time correlator Cq{t,t) = (l/N) {Sq{t)S^{t)) follows as 



Cq(i, t) = Cq(0, 0)Rl{t, 0) + 2T f dt' Rlit, t') (2.9) 

= + 2T f dt' 9jf)^-2^(t-t') (2.10) 

9{i) Jo 9it) 
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and we note for later the identity 



dtC^it, t)^2T- [2u;+^\ Cq{t, t) (2.11) 



git) 

The two-time correlator Cq{t,t„) = (l/N) {Sq{t)S^{t^)) can be deduced from the 
analogue of (|2.8|l for initial time 

S'q(t) = i?q(t, iw)5'q(<w) + / dt' i?q(i, t')^^{t') (2.12) 



as 

Cq{t,tvj) = i?q(<, iw)C'q(iw, iw) (2-13) 

The position-dependent correlation and response functions Cij{t,t^) and Rij{t,t^) 
are then just the inverse Fourier transforms of Cq(i,tw) and Rq(t,t^), respectively, 
with q conjugate to rj — r^. 

2.1. The function g{t) 

The calculations outlined above show that the Gaussian dynamics is fully specified 
once the function g{t) is known. The latter can be found from the spherical constraint, 
which imposes J{dq)Cq{t,t) = 1. Here and below we abbreviate (dq) = dq/(27r)'', 
where the integrals runs over the first Brillouin zone of the hypercubic lattice, i.e. 
q G [—tTjTt]'^. Using H2.10|l this constraint gives an integral equation for g{t): 

-2u}t 



g{t) = j [dq) Cq(0, 0)e-2-* + 2T dt' f{t - t')g{t') (2.14) 
where 

fit) = / (dq) = [e-^'lomf « (S^t)-'^/^ (2.15) 



Here Iq denotes a modified Bessel function and the final expression gives the 
asymptotic behaviour for large t. In terms of Laplace transforms g{s) = 
dt exp{—st)g(t), eq. (|2.14() then has the solution 

1 ^q(0:0) 



g{s) = ^ / (dq) ' ^ (2.16) 

With the exception of Sec.|Hl we focus in this paper on random initial conditions, 
Cq(0,0) = 1, corresponding to a quench at time t = from equilibrium at infinite 
temperature. In this case the q- integral in the last equation is just /(s), so that 

g{s) = (2.17) 

1 - 2r/(s) 

The asymptotics of the corresponding g{t) are well known; see e.g. O |^. For T 
above the critical temperature Tc, which is given by 

= 2/(0) = / (dq) - (2.18) 



there is a pole in g{s) at s = 2zeq- Here Zeq is found from the condition 2Tf(2zeq) — 1 



or 



/ 



{dq) — - — = 1 (2.19) 



2eq + 
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The presence of this pole tells us that g{t) ~ exp(2zoqt) for long times, implying 
that the Lagrange multiplier z{t) approaches z^q for f ^ cx). Correspondingly, the 
condition (|2.19|) is just the spherical constraint at equilibrium, bearing in mind that 
C^{t,t) = T/{zcq + Lu) from H2.5|l . Because lu k, for small q ~ |q|, the phase space 
factor in the q-integrals is [dq) ^ dwLo'^l'^^^ for small q or lj. This shows that as 
given by (|2.18|l vanishes as d ^ 2 from above; consequently we will always restrict 
ourselves to dimensions d above this lower critical dimension. 

At criticality (T — Tc), Zeq vanishes, and g{t) therefore no longer grows 
exponentially; instead one finds [HI I21j 

It is this case, of a quench to the critical temperature, that we will concentrate on 
throughout most of this paper. This is because here the FDR has the most interesting 
behaviour. 

We note briefly that, in principle, J{dq) should be written as {^/N)J2qi with 
the sum running over all q whose components are integers in the range —L/2...— 
1,0,1,... L/2 — 1 (assuming L is even) multiplied by an overall factor 2tt/L; there are 
N such q. When considering continuous functions of q this sum can be replaced by 
the integral J{dq) , and this will almost always be the case in our analysis. Exceptions 
are situations with a nonzero magnetization, where the wavevector q = is special 
and has to be treated separately. This is relevant in equilibrium below Tc, which we 
discuss briefly in Sec. 15.31 and for non-equilibrium dynamics starting from magnetized 
initial states (SeclSJ. 



2.2. Long-time scaling of Cq 

It will be useful later to have a simplified long-time expression for Cq(tw, ^w) for the 
case of a critical quench. At zero wavevector one has 

Co(tw, tw) = (l + 2T /*" dt' g{t')] « ^ (2.21) 

5(tw) V Jo J 1-K 

where the last approximation is based on H2.2()|l and is valid for large tw. For nonzero 
q, on the other hand, 

Cq(tw,tw) = +2T /*"dt'5(t')e-2-(*-*')) « -(2.22) 

Sl^wj V Jo J ^ 

which is as expected since all nonzero Fourier modes eventually equilibrate. The 

— 1/2 

crossover between the two limits takes place when ujtw ~ 1, or g ~ iw ; physically 

1/2 

this represents the growth of the time-dependent correlation length as ^ • We 
therefore introduce the scaling variable w — ut^: 

Cq{t^,t^) ^ e'^-' + 2Tuj-^ dyg{y/u,)e-^i-'-y) 
Co(tw,tw) I + 2Tlli^^ dy g{y/uj) 

Now keep w constant and let t^ — > oo, i.e. — > 0. Then g{y/uj) ~ (y/u)^^ and the 
second terms dominate in denumerator and nominator to give 

= i^-^)fdy y-«e-2-(i-^) (2.24) 

L/0(tw,Ewj Jo 



Fluctuation-dissipation relations in the non- equilibrium spherical ferromagnet 



8 



Combining with then gives the desired long-time scahng form 

Cq(tw,tw) - -^cMw), Tc{w)^2w( dyy—e~^^^^'y'> (2.25) 
^ Jo 

For d > 4: {k — 0) this simphfies to Tci^) = 1 — e^^™. As the derivation shows, 
eqs. (I2.24f) and H2.25|l are vahd whenever ^ 1, even for w = C(l). The latter case 
corresponds to w — > cxd and gives !Fc{w) — 1, which is indeed consistent with H2.22I) . 

For quantities such as Cq(tw,iw) that depend only on a single time variable, 
what is meant by the long-time limit is unambiguous. For two-time quantities like 
Cq{t,tv,) we use the following terminology: the long-time limit refers to the regime 
t :$> 1 and ^ 1 but without any restriction on t — t^, which in particular is 
allowed to be short, i.e. of 0(1). The aging regime indicates more specifically the 
limit tyj oo aX fixed ratio x ~ t/t^ > 1^ which implies that also t — t^ is large, of 
0{t^). Occasionally we specialize further to the regime of well-separated times, which 
corresponds tot:^t^:$> 1, i.e. the asymptotic behaviour of the aging hmit for x 1. 

To illustrate the difference, consider which wavevectors dominate the integral 
/ (dq) Cq{t, tw)- In the long-time limit at equal times t — t^, the scaling Cq(iw, ^w) ~ 
1/uj for uj ^ 1/iw combined with [dq) ~ du}Uj'^/^~-^ for small uj shows that the integral 
is divergent at the upper end of the frequency regime uj — 0{t~^) for all d > 2; in 
other words, it is always dominated by values of w (and therefore q) of 0{1). This 
remains true for two-time correlations, as long as t — tw = 0{1). In the aging limit, 
however, we have t — t^j — 0{t^) 3> 1 and the exponential factor from i?q in (|2.13|l 
then ensures that only values oi uo < (i — tw)^^ — 0{t~^) have to be considered in the 
integral. 



3. Fluctuation-dissipation relations for finite-range observables 

In this section we consider FD relations for observables that probe correlations over 
a lengthscale that can be much larger than the lattice spacing, but remains much 
smaller than the system size. The latter can then be taken to infinity independently, 
so that the 0(A^~^/^)-fluctuations of the Lagrange multiplier z become irrelevant. We 
begin by considering briefly spin observables, and then discuss bond observables in 
some more detail. 



3.1. Spin observables 



Since all observables that are linear in the spins can be written as superpositions of 
the Fourier modes S'q, the basic ingredient for understanding the FD behaviour is the 
FDR for the latter. Using (|2.11|l . this follows after a couple of lines as (C = dt^C) 

-I -1 



Xq{t, t^ 



TRq{t, tv 
C*q (t , tw 



= T 



2T- 



25 (iw) 



Cci (tw ; ^ W ) 



(3.1) 



This is independent of the later time t, a feature that is commonly observed in simple 
non-equilibrium models 2 . 

The fluctuating magnetization is simply Sq/N, so setting q = in H3.1|l gives 
directly the FDR for the magnetization 



Xo{t, tv 



T 



2T 



2g^t^ 



1 + 2T 



dt'g{t') 



(3.2) 
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As increases this converges on an 0(1) timescale to the limit-FDR 

yoo ^ T ^ = / 1/2 (d > 4) . . 

2T+{k/2)[2T/{1-k)] 2-k \ 1 - 2/d (d < 4) ^''■^> 

which is identical to the value obtained from the local magnetization as one would 
expect on general grounds. Without working out the susceptibility explicitly, it is 
clear from the t-independence of Xo{t,t^) and its fast convergence to X°° that the 
limiting FD plot is a straight line. Both of these observations are exactly as in the 
Ising model in d = 1 ^U]. Simulations have shown that also in the d = 2 Ising 
case the local-global correspondence holds for spin observables; the limiting FD plot 
is numerically indistinguishable from a straight line, though renormalization group 
arguments suggest that it should deviate slightly [THl [T71 ITT] . 

We should clarify that the Gaussian theory above applies directly not to the FDR 
for Sq but to the one for with q <^ iw"^^^ but q L~^, where L — N^/''- is the 
linear system size. The corresponding physical observable is a "block" magnetization, 
i.e. the average of the spins within a block of size £ ~ 1/g much larger than the time- 

1/2 

dependent correlation length ~ t„ but still small compared to the overall system size. 
For q = 0, i.e. £ = L, one would in principle need to account for the non-Gaussian 
fluctuations. However, it turns out that these are negligible as long as the system is 
not magnetized on average (see Sec.|Hl), so that the above results remain correct even 
for the global magnetization itself. 

More generally, the FDR for any finite-range spin observable can be expressed 
as a superposition of those for the Fourier modes; this can be seen by arguments 
parallelling those in the d = 1 Ising case ^Hl- As there, one can then show that the 
asymptotic FDR that is approached for well-separated times i ^ tw ^ 1 is dominated 
by the contribution from q = 0, and hence identical to X°° calculated above At 
equal times, on the other hand, equilibrated modes with q = 0{1) dominate and give 
X — 1. The crossover between these two regimes takes place when t — t^j — 0{t^) 
and follows (by superposition) from the corresponding crossover at g = 0{t^^) in the 
Fourier mode FDRs. From H2.25|l and H3.1(l the latter can be expressed as 

^q(i,iw) =^x(c^tw), T^\w) = 2-{2w-k) [ dyy-^e-^^^^-y^ (3.4) 

Jo 

in the long-time limit, providing the expected interpolation between X — X°° = 
l/[2 + k/(1 - k)] for w ^ and X ^ 1 ior w ^ oo. 



3.2. Bond observables 

Next we consider bond energy observables, ^{Si ~ SjY ■, where i and j are n.n. sites. 
Since all variables are Gaussian, the connected correlations follow by Wick's theorem. 
For the correlation of bond energies one gets 

Cy,M(i,tw) = 2^ ([5,;(t) - 5,(i)][5fc(tw) - Si{t^)\f = i [Cfc - Ca - + C^f (3.5) 

where time arguments have been left implicit. For the local case (i,j) = ik,l), this 
simplifies to 

Q,,., = 2[1 - a,]' (3.6) 

which tends to a nonzero constant for t — t^ ^ oo since Cy then approaches its 
equilibrium value, which is < 1. 
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Next we turn to the response function. In general, if one perturbs the Hamihonian 
H by —hB6(t ~ t^), then the equation of motion for Si acquires an extra term 
h[dB/dSi)5{t — tyj). So the perturbation in Si is 

f) H 

SS.,{t) ^ hyR,,it,t^) — {t„) (3.7) 



j 



dS, 



Thus the perturbation of an observable A is 
6A{t) = hJ2-g^{t)Rvit,t^ 

giving the response function |13j 



dB 



(^w) 



\ , ^ / dA dB 
RAB[t,t„) = y Rtj[t,t„){ TrTritjjTTritwj 



dSi ' ' dSi 



(3.8) 



(3.9) 



For A = 1(5-,, 



Sjf, B = i(S'fc - Sif this yields 

Rij,kliiy'tvj) — [Rik — Ril — Rjk + Rjl][Cik — Cu — Cjk + Cjl] (3.10) 

We now analyse the scaling of correlation, response and the resulting FDR. In 
terms of Cq{t, iw), the bond correlation H3.5(l is 

1 



Cij,kl {t, tvj) — T: 



(dq) Cq 



,jq-(ri-rfc) _ jq-(ri-r;) _ iq-(rj-rfc) , jq-(rj-r;; 



(3.11) 



We can take out a factor exp(jq • Ar) from all the exponentials, where Ar = 
i(rj + Tj) + i(rfe — r;) is the distance vector between the bond midpoints. In the 
remaining exponentials, q is multiplied by vectors with lengths of order unity. 

Now assume t — t^^l. As explained above, integrals of two-time quantities over 
q are then dominated by the small-g regime, « w < (i — tw)^^- We can therefore 
Taylor expand in q and get, using the equivalence of the d lattice directions, 

n 2 



Cij^ki — 7: 



d-\r^ 



Tj) • (rfe -r;) J {dq)Cqq 



2giqAr 



(3.12) 



(3.13) 



Similarly one finds for the response 

R^J,kl = [d-\v,^Vj)-{vk~ri)f j{dq)R^q^e'^'^'' j {dq) C^q^e'^'^'' 

For the local bond-bond correlation and response one sets Ar ~ and has (r^ — r^) 
(r^ — Tj) = 1, which gives for the FDR 

J{dq)TR^q^ _ J{dq)R^q^ 



-'^bondC*' *w) 



J{dq)C^q^ J{dq)X^'R^q^ 



(3.14) 



So 1/X^°^j(t, tw) can be thought of as an average of X^'^{t,tvj) over q, with the 
weight Rq{t,tw)q^^ ■ The factor i?q ensures that significant contributions come only 



from wavevectors q up to length g ~ {t — t^, 



-1/2 



, i.e. up to uj{t — tw) ~ 1- Thus, 



when t~tvj ^ iw, the result is dominated by the regime ujt^ ^ 1, where Xq ~ 1. For 
t — tv, ^ iw, meanwhile, one only gets contributions from ujt„ <^ 1, where Xq = X°° . 
So the FDR 13.14|l for local bond observables is a scaling function interpolating 
between 1 and X°°, with the same X°° as for the magnetization. Explicitly one 
has, using (|2.()l) and changing integration variable from q to w = ujt^, 



"v^loc 
^bon. 



d(^' ^w) 



(3.15) 
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with X — t/tyj and Tx the scahng form H3.4|) of Xq. To find the shape of the FD plot, 
recall that the equal-time value of the local-bond correlation (|3.t)|l is a constant in the 
long-time limit. For t - 
as 



^ 1, on the other hand, eq. ^A.VZ\ shows that C\°^^^ scales 



9{U 



git) 
git) 



dujLu'^/^-^-Tc{c^t^)e-^^'-'-'>L 



1 2 



t: 



dww'"^-'Tciw)e-^^'-'-'^"- 



(3.16) 



(3.17) 



Since Tciw) — > 1 for w oo, the w-integral would be divergent without the 
exponential cutoff and scales as [it — t^)/t^]~'^/'^ for t — t^ <^ t^, so that C^°nd(^' tw) ~ 
{t — t^)~'^ in this regime. The regime t — > where -'^bond(*' ^w) 7^ 1 is therefore 
compressed into the region where C]^^^-^^ is of order t~'^, so that the long-time limit of 
the FD plot is a straight line with equilibrium slope. Qualitatively one thus has the 
same behaviour as for local bond observables in the Ising model |10| . 

Next consider long-range bond observables, where we sum (ij) and (kl) over all 
bonds. The same proviso as above for the magnetization applies here, i.e. by applying 
the Gaussian theory we are effectively considering the bond energies averaged over 
a block that is large but has to remain nonetheless small compared to the system 
size. One can show that the resulting equal-time correlation again approaches a 
constant value for t^ — > cxi. (This follows because for large Ar, one can use the 
small g-expansion H3.12|l even for equal times. From Cq{t^,t^) w T/q"^ one gets 
C(Ar;t^,tw) ~ |Ar|2-'' for large Ar and so J{dq) Cqq^e'^i-^" - V^jArp-'' ~ |Ar|"''. 
The square |Ar|~^'^ then yields a convergent sum over Ar.) So we focus directly on 
the regime t — t^, ^ 1, where the expansion (|3.12|) is again valid. Keeping the bond 
(ij) fixed, the scalar product (r^ — Vj) ■ (r^ — r;) means that only bonds {kl) parallel 
to (ij) contribute, so that the sum over (kl) becomes a sum over Ar, running over 
all lattice vectors. (For non-parallel bonds, Ar could also assume other values not 
corresponding to lattice vectors.) The sum over (ij) then just gives an overall factor 
of Nd. Normalizing by N, the block bond correlation function is 

-I 2 



x^block/i , \ 
'-'bond ^vj) 



--y 

2dN ^ 

Ar 



idq) Cqq e 



2 iq-Ar 



1 

2d 



and similar arguments give for the (normalized) response 



oblock 



bond (''I ^ 



(dq) RqCqq 



{dq)C'y (3.18) 



(3.19) 



so that the FDR becomes 



block 



(i,tw) 



J{dq)TR^Cy _ J{dq)R^Cy 



Jidq)X^ 



(3.20) 



Again, this is the inverse of a weighted average of ^, now with weight RqCqq^. 
The same arguments as for (|3.14|l then show that X^'^°^'^{t,tyj) scales with x = tjt^ 
and interpolates between X = 1 for x = 1 and X = for x — + 00. The value of 
the correlator H3.18|l decays from 0(1) at t = to 0(t^'^^'^) at the point x — \ k, \ 
where aging effects appear. While this is larger than for the local bond observables, 
it still decreases to zero for t^ (x>,so that the limiting FD plot is again of pseudo- 
equilibrium form. This is different to the case of the Ising model, where the global 
bond observables give nontrivial limiting FD plots |10|. 
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In more detail, the scaling of the block bond correlator ()3.18|l in the aging regime 
t - > 1 is 

«(i,tw) - / dcou''/'-'^TUcot^)e~'-^'~'-W (3.21) 

9{t) J uj^ 

^ (^Ly t-^/^ jdww^/^-^TUw)e-^'^^-^'^^ (3.22) 

The integral scales as {x - l)-"^/^ for x w 1, so C^o°nd (t,iw) ^ {t - t^y^/"^ there. 
For x ^ 1, on the other hand, the integral becomes ~ ^-{d+i)/2^ ^^biock^^^ ^ 
^_«+2^K-(d+4)/2^ Explicitly, in this t > t^ regime, C^^°'f - <^t-('^+4)/2 for d > 4, and 
Cboiid^ ^ ti/'^t^'^ for d < 4. The response function scales in the same way as dt^C^^^j^^, 
because X is everywhere of order unity. 

3.3. Product observables 

Instead of the bond observables ■^{Si — Sj)^ we could consider the spin products 
A — SiSj, B = SkSi. The correlations are then 

C„-,fci(i,iw) - {S^{t)Sj{t)Skit^)Si{t„)) - {S,{t)Sj{t)) {Sk{t^)Si{t„)) (3.23) 

= akit,t^)Cjiit,t^) + Ca{t,t^)Cjk{t,t^) (3.24) 

The local equal-time correlation function Cij,ij{t^,t^) thus approaches 2 for oo. 
The corresponding response function is 

Rij,ki{t,tw) = RikCji + RilCjk + RjkCil + RjiCik (3.25) 

In the ZocaZ case, one can replace all functions by local ones in the aging regime - there 
are no cancellations leading to extra factors of as was the case for bond observables, 
compare e.g. (|3.11|) and H3.12|l - so that the FDR 

^prodV''' ''w^ 4C'C - C 

becomes identical to the one for the spin autocorrelation and response. In particular, 
one again gets the same X°°. 

For the global (block) case, we can write 

Cpro?(^)^w)= ^ Cij^kl{t,tv,) (3.27) 

fe),(fcO 

= [Crk{t,t„)Cji{t,t„) + Cu{t,t^)Cjk{t,t„)] (3.28) 

ijkl 

= ^IZ"»j"fei'^^fc(*'*w)Qi(i,iw) ^^J {dq)n\Cl{t,t^) (3.29) 

ijkl 

where riij = 1 if i and j are nearest neighbours and otherwise, and riq is its Fourier 
transform. For the response one has similarly 

^prof ^w) = j (dq) nlR^it, t^)Cq{t, iw) (3.30) 

In the aging regime, where t — tw S> 1, the integrals are dominated by small q, where 
riq can be approximated by the constant no = 2(i. This cancels from the FDR, which 
becomes 



yloc (4- 4- \ _ ^^cJ^^ii'^ii _ ^cJ^Hi /„ „„N 



x^block/. . \ _ /(rfg) TcRqCq 
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This is the inverse of the average of -'^"q ^ with weight RqCq. Again, this is a scahng 
function of t/t^ interpolating between 1 and the same as for spin observables. 

The scahng of the block product correlation function (|3.29f) itself is a little more 
complicated than for the bond observables and depends on dimensionality. Focussing 
again on t - tw > 1 one has CpJ°Jj'^(t, tw) w 2d'^ J (dq) C^{t, t^). The integral defines 
the function CC{t, t^i) discussed in Sec. for d > A and Sec. [3 for d < 4. In the 
former case, one has from H6.8H6.10| that CC{t,tw) = CCcq{t — t^)Tcc{t/tw) where 
CCoq(i — tw) ^ {t — tw)^^~'^^^^ asymptotically; this equilibrium contribution governs 
the behaviour of Cp^°J^'^(i, t^) for t—t^ <C t^. Where aging effects appear {t — t^ ~ t^), 

^prod' ~ tw~'^^^^ and so one gets a limiting pseudo-equilibrium FD plot. In the regime 
of well-separated times x 1, the scaling function Tcc{x) decays as x""^ so that 
Cprod'(^'*w) ~ tlt-'^/'^. These scalings, though not the overall magnitude of C^pl°f, 
are the same as for the energy correlation function Ce in (|6.14|) below: both functions 
are proportional to CC{t,t^) in the aging regime. 

In the opposite case d < 4, the equal-time value CC{t, t) (and there- 
fore Cpl°'^{t,t)) diverges as see ()7.4|l . The normalized correlator 
CC{t,tv,)/CC{t,t) is a scaling function G{x) of a; = t/tw, implying that the nor- 
malized FD plot will approach a nontrivial limit form, with asymptotic slope 
as shown above. Quantitatively, because G{x) ~ x~'^l'^ for a; ^ 1, one has 
0'(^'*w) - t'^^-'^y^t^/tr/' ~ tTt'"' for t » 



4. Correlation and response for global observables 

We now ask what happens if we go from block observables to truly global ones, which 
reflect properties averaged over the entire system; the total energy is an important 
example. We anticipate that here non-Gaussian fluctuations are important. Indeed, 
the results above show that this must be case. Otherwise we could directly extend the 
Gaussian theory results from block to global observables, with no change to correlation 
and response functions. The global bond observable is just the energy. Using the 
spherical constraint, this can be written as 

^ = E ^(^^ - ^j)' = ^ ~ E ^^^^ (4.1) 

(ij) {ij) 

and so is identical to the global spin product observable, up to a trivial additive 
constant and sign. So the global bond and product observables must have identical 
correlation and response functions; but we saw above that this requirement is not 
satisfied by the Gaussian theory. Thus, non-Gaussian fluctuations are essential to get 
correct results for global observables. 

Physically, the origin of the distinction between block observables and global 
ones is the effective infinite-range interaction induced by the spherical constraint. In a 
model with short-range interactions, block observables will show identical behaviour to 
global ones whenever the block size I is larger than any correlation length in the system, 
whether or not £ L: the behaviour of any large subsystem is equivalent to that of 
the system as a whole. In the spherical model, the infinite-range interaction breaks 
this connection, and global correlation and response functions cannot be deduced from 
those for block observables. 
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4.I. Non-Gaussian fluctuations 

To make progress, we need to return to the original equation of motion 1)2. 2|l . This 
can be written as 

dtS, = - ^ n,,Sj + 6 - {z{t) + N-^/^Az)S, (4.2) 

j 

where the notation emphasizes that the fluctuating contribution to the Lagrange 
multiplier is of 0{N^^^^). The latter induces non-Gaussian fluctuations in the Si 
of the same order. This shows quantitatively why the Gaussian theory works for 
block observables: as long as one considers correlations of a number of spins that is 
< N, fluctuations of ©(TV^i/^) 

can be neglected. For global observables, on the other 
hand, we require the correlations of all N spins and the Gaussian approximation then 
becomes invalid. 

To account systematically for non-Gaussian effects we represent the spins as 
Si — Si -\- N~^^^ri, where Si gives the limiting result for N — > oo, which has purely 
Gaussian statistics, and N~^^^ri is a leading-order fluctuation correction which will 
be non-Gaussian. Inserting this decomposition into 1)4.2(1 and collecting terms of 0(1) 
and ©(TV"^/^) gives dtSi — — flijSj — z(t)si-{-S,i as expected; to lighten the notation we 
use the summation convention for repeated indices from now on. For the non-Gaussian 
corrections one gets the equation of motion 

dtri = -^ijrj - z{t)ri - Az Si (4.3) 

with solution 

r,{t)^R,j{t,0)rj{0)~ f dt' R,,{t,t')sj{t')Az{t') (4.4) 

'0 



The properties of Az{t') can now be determined from the requirement that, due to 
the spherical constraint, N^^ J2i ^ii^) = 1 at all times. Inserting Si — Si + N^^/^ri 
and expanding to the leading order in iV^^/^ gives the condition 

1 s,it)r,{t) = -\n-''^ Y.^sU^) - 1) ^ -^^W (4.5) 

i i 

where the last equality defines A(i), a fluctuating quantity of 0(1) that describes the 
(normalized) fluctuations of the squared length of the Gaussian spin vector Si. At 
t = Q, the condition H4.5() is solved to leading order by setting ^^(0) = — iA(0)si(0), 
since (1/A^) (0) = 1 + 0(iV"^/^). With this assignment, and setting a{t) = 
2Az{t) + A{Q)5{t), eq. 631 reads 

n{t) = J dt' R,j{t,t')s,{t')ait') (4.6) 

We have left the integral limits unspecified here: the factor Rij automatically enforces 
t' < t, and we use the convention a{t') = for t' < 0. The spherical constraint 
condition (|4.5I) then becomes 

dt' ^s.,{t)R,,{t,t')s,{t')a{t') = A{t) (4.7) 

Now, up to fluctuations of 0{N^^^^) which are negligible to leading order (even if 
they are correlated with a(t')), we can replace {l/N)si{t)Rij{t, t')sj{t') by its average 

K{t,t') = 1 {s,it)R,,{t,t')sjit')) = 1-R^^(t,t')a,{t,t') = I (dq) R^it,t')C^{t,t'){A.8) 
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If we then define the inverse operator, L, of K via 

dt' K{t,t')L{t',t„)^S{t-t^) (4.9) 
for tw > 0, then the solution to H4.7() is 

a{t) = / dt' L{t,t')A{t') (4.10) 



where for consistency we adopt the convention A(t') = for t' < 0. With (|4.()|) we then 
get an explicit expression for the non-Gaussian ©(A^^^/^^.^orrections to the spins, 

nit) I dt'dt" R,,{t,t')s,{t')L{t' ,t")A(t") (4.11) 



2 . 

in terms of the properties of the uncorrected Gaussian spins Si . 
4-2. The functions K and L 

Before proceeding, we analyse the properties of K and L. From 14. 811 . K{t, t') vanishes 
for t < t' while its limit value for t t'+ is {l/N)d^jC,j(t,t) = {l/N)CH(t,t) = 1. 
Inserting (ESJ, (|TTT1) and itTT^ into l|TH|l . one also finds that the equal-time slope has 
the simple value dfK (t, t')lt=t'+ = ^T. From these properties and the definition l|4.9|l 
it follows that 

L{t,t') ^S\t-t') + L''^\t,t'), L'^^Ht,t') ^2Td{t^t')- L'^^\t,t') (4.12) 

where L'-^\t,t') vanishes for t < t' and jumps to a finite value at t ~ otherwise 
it is smooth and, as we will later see, positive. The structure of (I4.12|l can be easily 
verified e.g. for the limit of equilibrium at high temperature T, where z^q = T and all 
ui can be neglected compared to Zeq. One then has K{t,t') — exp(— 2T(i — t')) and 
the inverse (|4.9|l can be calculated by Laplace transform. Since the Laplace transform 
of K{t, t') is K (s) = l/(s + 2T) this gives L{s) ^ s + 2T, which corresponds to 
with ^ 0. 

We next determine the long-time forms of K and L'^) foj- quenches to criticality. 
In both cases it is useful to factor out the equilibrium contribution. For K this is, 
from ()4.8|) and using H2.6|l and (|2.13|) . 

K,q{t - tw) - / (dq) e-2-(*-*w) (4.13) 
J 1^ 

Apart from the factor of 2 in the time argument, this is just the (critical) equilibrium 
spin-spin autocorrelation function. One can also write K{t) — 2Tcj^ dt' f{t') 
from (|2.15l) and this shows that Kcqit) ~ i(2-'')/2 fQj- large time differences. The 
ratio K{t,tyf)/ Keq{t — t„) will show deviations from 1 when aging effects appear, i.e. 
when t — t^ ~ t^. The form of these deviations can be worked out by using the scaling 
form H2.25|) of Cq(t.jv, t^) and recalling that only the small g-regime contributes, where 
{dq) ^ dio uj'^/'^^^ . Changing integration variable to w — Lut^, gives 

K{t, iw) = K,q{t ~ t^)TK{t/t^) (4.14) 

, fdu;w;('^-4)/2e-2(--i)-^c:(w) . ^ 

^-(-) =-^ ' J,^^(.-4)/2,-2(.-r); (4-15) 

where the first factor in 14.15|) arises from the two factors of [g(iw) / git)]^^"^ contributed 
by Rq{t,t^) and Cq(i,tw), respectively. By construction, J-'k{x) should approach 1 
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for x — !■ 1; indeed, in this limit the lu-integrals are dominated by large values of 
w ~ — 1), for which J-c = 1. The decay for large x follows from J-'c{w) ~ w for 
small w as J-k{x) ~ x'^^^ . Explicitly, one finds by using (|2.25|l and carrying out the 
ly-integrals that 

^k{x) = ^—^{x - l)(d-2)/2^. dyy'^ix - y)-''/^ (4.16) 
2 Jo 
For d > 4, where k = 0, this gives 

J'k{x) = 1- (d>4) (4.17) 

while for d < 4 the required indefinite integral is [{d— 2)/2] J dy y'''^~^^^^{x — y)~'^^^ = 
—x~^{x/y — l)(2-'^)/2 and one gets simply 

Tk{x) = a;(2-'i)/2 (rf < 4) (4.18) 
Next we determine L^^^. Combining (|4.9|) and (|4.12|l . the defining equation is 

dt' K{t, t')L^^^ {t', iw) = 2TK{t, t^) - dt^K{t, t^) (4.19) 
Again it makes sense to extract the equilibrium part of L^'^\ This is defined by 

dt' K,^{t ~ t')L^J'^ it' ~ t„) = 2ri^eq(Ai) + K^i^t) (4.20) 
where lS.t ~t — t^. Solving by Laplace transform gives 



^(.)(^)^2Ti:e,(.) + (^i^c,(.)-l)^^^^,^__^ 

where from (|4.13ll . at criticality, 



(4.21) 



i^e,(.)=r./(rf,)-^ (4.22) 

The leading small-s behaviour of this is -ftreq(O) — K^qi^s) ~ s(''-4)/2 for d > 4 (plus, 
for d > 6, additional analytic terms of integer order in s which are irrelevant for us). 
For d < 4, on the other hand, i^oq(s) ^ s(''~4)/2 divergent for s — > 0. Inserting 
these scalings into (|4.21ll and inverting the Laplace transform gives for the asymptotic 

(2) 

behaviour of Lcq 

^eq {t) ^ I ^(,_6)/2 ^ 4) (4.23) 

f2l 

It will be important below that, for d > 4, Kcq{t) and Lcq {t) both decay 
asymptotically as t(2-(i)/2_ r^^ie ratio between them can be worked out from (|4.21|) . 

by expanding for small s as l/-ft:cq(s) « l/[i^oq(0) - cs^"^-^^/^] = l/i^oq(0) + 

cs(''-4)/2^^2^('Q') .ftrhere c is some constant; comparing with i^eq(s) ~ -R^eq(O) — 

^g{d-4)/2 gjygg 

L?q)(t)=ifeq(i)/i^e'q(0) (4.24) 

for large time differences t. 

The integral of (t) over all times follows from H4.21|) as 

ig.(o)=jr".,4;.(<)={ g-*-<''*^'^=i'-'/-f''''''<^=/'^" (425) 
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Using the fact that J{dq) {T^/uj) 1, one has /(dq) {T^/ujf > 1 so that Lf^{{)) is 
positive independently of d. This is consistent with the intuition that, with the sign 
as chosen in (I4.12|) . the function is positive. 

With the equihbrium part of L^"^^ determined we make a long-time scaling ansatz 
for L(2), 

L(2)(t,t^) = (i - t^)TL(t/t^) (4.26) 
so that H4.19|l becomes 

dt'K,^{t - t')4q^(i' - t^)^K{t/t')TL{t' /t^) = 

- 2T,K,^{M)Tk{x) + K'^^{M)Tk{x) + ^K,^{^t)T'K{x) (4.27) 

where At — t — t^ and x = t/ty„ as before. We now take the aging limit of large t^ 
with At = 0{t^) to determine !Fl- The second and third terms on the r.h.s. are then 
smaller by factors of order \/t^ than the first, and can be neglected to leading order. 
The second term on the r.h.s. of l|4.20|l is likewise subdominant, and this can be used 
to rewrite the dominant first term on the r.h.s. of (|4.27|) . giving 

j,^ ^ jdt'K^S - t')Lil\t' - t^)TK{t/t')TL{t' /t^) 28) 
J dt' K,^{t - t')Li^^ {f - t^) 

(2) 

We consider first d > 4. Then both the functions Kcq{At) and Lcq (At) have finite 

^ ^ (2) 

integrals Kcq{0) and Lcq (0), respectively, over At = . . .00. In the aging limit, the 

(2) 

factors Kcq{t—t') and Leq (t'—t^) therefore act to concentrate the mass of the integrals 
appearing m around t' = t„ and t' = t. This can be seen more formally by 

changing to y = t' /t^j as the integration variable and taking t^ large. Then the factors 
Kcq{tyj{x — y)) and Liq(tvj{y — 1)) produce singularities ~ (a; — ?;)(2-'*)/2 for j/ ^ x 
and ~ (y — 1)(2^'')/-^ for y 1, respectively, and because these are non-integrable 
they dominate the integral for tw — > 00. All other factors in the integrals are slowly 
varying near the relevant endpoints and can be replaced by their value there. In the 
aging limit we can therefore write H4.28(l as 

^ Koqmi'^{At)TK{imx) + K,q{At)W^iO)TK{x)TL{l) 29) 

i^eq(O)i^qHAt) + Xeq(At)L('q^(0) 

Eq. (|4.24l) tells us that the At-dependent factors cancel, giving together with (|4.25|l 
and.FK(l)=^L(l) = l 

^ , , K-q\0)TL{x) + Li^^iO)TK{x) J^L{x) + [2T,k,qiO)-l]J^Kix) ,^ 
■Fk[x) — — ^ — ; ^r7T\ = (4.30) 

ifcqnO)+4q^(0) mK,q{0) 

In (i > 4, where i^cq(O) is finite, we therefore have the simple result that the scaling 
functions for K and L'^^^ are identical, 

Tl{x)^Tk{x) (4.31) 

But in the limit d — > 4 from above, K^qiS^) diverges and H4.30|l gives no information 
about Tl- For d < 4 a different approach is therefore needed to determine Tl- One 
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subtracts from (|4.27|) the first and second terms on its r.h.s., using (|4.2U|) to rewrite 
them as an integral and changing integration variable from t' to y — t' jt^. This gives 



1 



-K,4^M)T'k(x) (4.32) 



dyK,q{t^{x - y))Lf^{t^{y - 1)) [Tk {x/y)TL (y) - ^A-(a:)] 

_ t 

For y ^ X, Kcq{tw{x — y)) contributes a singularity ^ (x — y)^^"''^/^ which is integrable 
in d < 4. For y — s- 1, the terms in square brackets vanish as ~ y — 1 since J^Liy) is 
smooth at y = 1 as we will see below, in the sense that .?i(l) is finite. These terms 
combine with the ~ (^ _ l)(^-6)/2 fj-Q^-^ ^^^^ integrable (y — l)'''^'*)/^. The 

(2) 

contributions from the short time behaviour of Kcq and Lcq are therefore unimportant 
in the aging limit and we can replace these functions by their power-law asymptotes. 
Up to overall d-dependent numerical factors the condition H4.32|l then becomes 

t-' j'^ dy {x - y)^'-'^/\y l)''^-^)/^ [:Fk {x/y) {y) - J^k^x)] = 

= ±At(2-'^)/2^^(^) (4.33) 

In the aging limit At scales as t^, and so the l.h.s. of this equation (^ t~^) becomes 
large compared to the r.h.s. {'^ tw''^^) unless the y-integral vanishes. The required 
condition for J-^ is therefore 

dy{x - 2/)(2-'^)/2(2/ - 1)^'-'^/' {Tk {x/y)TL (y) - Tk{x)] = (4.34) 



This is in principle an integral equation for Tl- Fortunately, however, the solution is 
the naive extension of H4.31|l to c? < 4: with Tk{x) = x'^^"''-'/^ from H4.18|l one sees 
that for !Fl{x) — J-k{x) — x^^^'^^'l'^ the square bracket in 14.34|l vanishes identically. 
The identity H4.31|l therefore holds both for d < 4 and for d > 4. 

In summary, we have determined long-time scaling forms for K and L for quenches 
to criticality. For K, the result is (|4.14|l with (|4.13|l and (|4. 1714. 18(1 : for L, we 
have with l|4. 2114. 23(1 and (|OT|l . Combined with this fully determines 

the leading non-Gaussian corrections to the spherical model dynamics (at long times, 
and after a quench to criticality from a random initial state). 

5. General expressions for energy correlation and response 

In this section we derive general expressions for the two-time correlation and response 
functions of the energy, taking into account non-Gaussian fluctuations. The results 
will be valid for arbitrary quenches since we will leave K and L unspecified. 

5.1. Energy correlation junction 

We can write the energy 1(2.1(1 as H — ^SiflijSj. Inserting Si — Si + N^^/^ri, the 
energy correlation function (normalized by N) is to leading order 

CEit,t^) ^ ^ (^(^S^n^JSJ +2N-^/^S^n^JrJ^ ^ (^SkflklSl + 2N-^^hkflkiri^ ^ ^ (5.1) 

Using ((4.11() , all quantities involved can be expressed in terms of the Gaussian variables 
Si so that the average can be performed using Wick's theorem, i.e. by taking products 
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Figure 1. Illustration of Wick pairings for . (a) The six Gaussian spins that 
need to be paired in 15.81 are indicated by circles with site labels. Spins arising 
from the expansion of H{t) and H{t^) are to the left and right of the vertical 
dashed line, respectively; any pairing contributing to the connected correlator 
must have links across this line. Dotted lines connect spins that are already 
coupled in space: Si and Sm via the factor QijRj^, and Sfc and s; via Qkl- The 
vertical lines between s„ and Sn indicate that pairings which couple these spins 
are not allowed, (b) The solid lines show the only Wick pairing that contributes 
to leading order in 1/A'^: it gives two independent groups of spins, (c) only gives 
one group and so is subleading. (d) has two groups but is excluded from the 
connected correlator because there are no pairs across the dashed line. 



of all possible pairings. We use the prime on the average to indicate the connected 
correlation function. This just means that in the Wick expansion all terms not 
containing any pairings of a variable at t with one at have to be discarded, 
since these terms give the disconnected contribution )• Multiplying 

out 1)5. l|l one obtains four contributions. The first one is 

1 2 

4cj^'^ = —n^jfiki (s,(t)sj(i)sfc(tw)s/(iw))' = ^^^y■^^fc^Qfe(^,^w)c,^(^,^w) (5.2) 

To eliminate one of the factors of fl, note from (|2.6|l that {dt + z{t))Rq(t,t^) — 
—ujRqit,t^) for t > t„. In real space, this reads 

{dt + z{t)) R,k{t,t^) = -n,,Rjk{t,t„) ^ -R,j{t,t^)n,k (5.3) 

and from H2.13|l an exactly analogous relation holds for Cij{t,tv,)- Thus 

^C^E^ = - ^[{dt + z{t))C^k{t,t^)]^kiCa{t.t^) (5.4) 

= - ^{dt + 2z{t))C,k{t,t„)^uCa{t,t^) = ~{dt + '2z{t))C^C{t,t^) (5.5) 

The last equality defines CilC, which is just the normalized trace of the product 
of the matrices Cik{t,t^), Qki and Cu = Cu; in Fourier space, CVLC{t,t^) = 
j{dq)ujCl{t,t^). 
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The second contribution to Ce is 

^4'^ = (s.Wr,(t)sfc(t„)sKtw)>' (5.6) 

= - ^ J dt'dt" n,jnkiR,m{t,t')L{t',t") 

X (s,(t)s„(t')iV"'/'A(t")sfc(tw)s;(tw))' (5.7) 
where we have inserted H4.11|l . Writing iV^^/^A(t") exphcitly this becomes 

4Cj,')= dt'dt" n,^nkiRjM')L{t\t") 

X (5,(i)s„.(t') [slit") ~ (slit"))] SkitMt^))' (5.8) 

We now need to perform the Wick expansion of the average. The subtraction sl — l^s^ 
means that ah terms which would pair s„ with s„ are excluded; Sk and Si also cannot 
be paired because we are considering the connected correlation. We can reduce the 
number of pairings further by considering that we need to get an overall result of 
C(l). The index j does not need to be considered further: after summing over j, 
is some translationally invariant function of the distance vector between spins 
Si and Sm- If the remaining indices i, k, I, m, n where unrestricted, then together with 
the 1/7V^ prefactor we would maximally get an 0{N^) result. Each of the factors 
^ijRjm and VLki couples two indices and so reduces the order of the result by 
Having already got two such couplings outside the average, we can only "afford" one 
extra coupling from the Wick pairings to get a contribution of 0(1). After some 
reflection one sees that this only leaves the pairing [im][fcn][Zn]: [im] introduces no 
further coupling beyond QijRjrm and [fcn][/n] gives only one further coupling beyond 
flki- Alternatively, we can think of this pairing as having the indices i, m and fc, I, n, n 
in two independent groups; each group gives a factor of N and this just cancels the 
1/iV^ prefactor. All other allowed Wick pairings give smaller terms, as illustrated 
graphically in Fig. njL,b. For example, [ifc][mn][ln] together with flijRjm and flki 
couples all indices into a single group and thus gives a term of only 0{1/N) (Fig.Q;). 
The pairing [m][mri][fcZ] would give two independent groups and thus an 0(1) term, 
but is excluded because k and I cannot be paired in the connected correlator (Fig. QJi) . 
Bearing in mind that our dominant pairing has a symmetry factor of 2 because the 
s„'s in [fcn][/ri] can be swapped we have thus finally 

4cf)= -2 J dt'dt" j^n,,Rjrn{t,t')c,ut,t'W,t") 

X ^nklCnk{t",t^)Cnl{t",t^) (5.9) 

dt'dt" [{dt + 2z{t))K{t, t') - 5{t - t')]L{t', t")CnC{t", t^) (5.10) 

In going to the last line we have exploited H5.3(l to eliminate ftij. Since t' is an 
integration variable, we have also been careful here to subtract off with the S{t — t') 
the spurious contribution which the dt applied to the step discontinuity in K(t,t') 
would otherwise give. The ^'-integration can now be carried out using H4.9|l and we 
get 

4cf ^ = {dt + 2z{t))cnc{t, iw) - / dt'L{t, t')cnc{t', t^) (5.11) 
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Figure 2. Wick pairings for . (a) represents the constraints for the possible 
pairings in I5.16i . (b) is the only pairing that contributes to leading order, forming 
three independent groups of spins. 



2z{t)cnc{t,t^)- dt'L^^\t,t')cnc{t',t^) (5.12) 



Note that we could use the same trick here as in (|5. 315. 5(1 to write CQ,C{t' ,t^) — 
— {\dti + z{t'))CC{t' ,t^), with CC{t,tyj) defined in the obvious manner as CC = 
N~^CikCki — J{dq)ujC^. However, this reduction from CQC to CC apphes only for 
t' > tvr, while for t' < t^ one would need to take a time derivative w.r.t. t^ instead of 
t'. This case distinction would make evaluation of H5.12|l awkward, so we retain CftC 
here and below. 

The third contribution to C^; is obtained by simply swapping the roles of t and 
in (|5.12|) and remembering that CflC is symmetric in its time arguments, 

4^"^ = J^^^i^hi {s^it)s^{t)sk{tMt^)y (5.13) 

= 2z{t^)cnc{t, t^) - I dt'L^^^t^, t')cnc{t, t') (5.14) 



The fourth and last contribution to Cb is, using again H4.11|l and writing out 
A(i") and ^{t'i) 

= -^^^J^ki {s^{t)rJ{t)sk{tMtw)y (5.15) 



dt' dt" dt'^dt'^ D,ijQkiRjmit, t')L{t', <")i?/„(t„, t'^)L{t'^,t'lj) 



1 

X (s,(i)s,„(t')s'(i")sfe(iw)s„(C)Sp(C))' (5.16) 
where it is understood that, because of subtractions which we have not written 
explicitly, pairings of So with itself and of Sp with itself are to be excluded. The 
only pairing that gives an overall C(l) contribution turns out to be [im] [fcn] [op] [op] 
and gives the required 3 independent index groups (i,m; A;,n; o, o,p,p) to cancel the 
prefactor; see Fig. |21 With the symmetry factor 2 for the possible swap of the 
[op] [op] pairings one gets 

4^"^ = 2 / dt'dt"dt'^dt'ij^n,^R,^{t,t')c,u-t,t') 



X -^^klRlnitw, t'^)Ckn{tw, t'^)L{t', t")L{t'^, C)^C'op(t", t'^) (5.17) 

dt' dt" dt'^dtl [{dt + 2z{t))K{t, t') - 5{t - t')]L{t\ t") 



2 , 

X [{dt„ + 2z{t^))K{t^, C) - 5{t^ - t'J]L{t'^, t'i)CC{t", C) (5.18) 
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Using H4.9|l one can carry out two of the time integrations to get 



4Cjj') = i y dt'dt'^ [2z{t)S{t - t') - {t, t')] 

X [2zit^)6{t^ ~ 4) - C)]CC(i', C) (5.19) 



5. ^. Energy response function 

To find the energy response, consider perturbing the Haniiltonian by a term —hd(t — 
tv,)H = —{h/2)5{t — t^)SiVLijSj, where h is the field conjugate to the energy. The 
equation of motion in the presence of the perturbation is therefore 

dtS., = -n.,jSj - {z{t) + hAz{t))S., + hS{t - t^)nijSj (5.20) 

where hAz is now the change in the Lagrange multiplier induced by the perturbation. 
The fluctuating component of z of 0{N~^^^) is in principle still present, but negligible 
compared to hAz for field strengths that are 0(iV°). Inserting a corresponding 
expansion of the spins, Si = Si + hri, gives for Si the unperturbed equation of motion 
and for the perturbed component 

dtri = -^ijrj - z{t)ri - Az{t)si + 5{t - t^)QijSj (5.21) 

The solution of this is 

r^{t) = R^k{t,t^)nkisi{t^) - I dt' R^k{t,t')Az{t')sk{t') (5.22) 



One now needs to determine Az. This can be done by noting that the normalized 
length of Si is 

j^T.Sf-j^T.^^ + s.r. + Oih^) (5.23) 

i i i 

The change to first order in h must vanish to preserve the spherical constraint, giving 
the condition (1/A^) J2i i^i'^i) = or, using (|5^^ . 

1 . .... /•*.., 1 



^-i^,fe(^,^w)^^fe^a^(^,^w) = J dt' —R,kit,t')akit,t')Az{t') (5.24) 

In the integrand one recognizes the definition H4.8|l of K, so that one can write the 
solution of this as 

Az{t)^ J dt' L{t,t')RnC{t',t^) (5.25) 

with obvious notation for RilC. 

Now we can find the change in the energy, 1/(2A^) (SiilijSj), which is given by 
(h/N) {ri^ijSj) to hnear order in h. Dividing by h and using H5.22|l then gives the 
energy response function 

RE{t,t^) = j^{r,n,jSj) (5.26) 



—Rik{t, t^)nki^ijCji{t, tyj) 

dt'dt" RnC{t, t')L{t', t")RnC{t", iw) (5.27) 
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One can eliminate one of the time integrals by using (|5.3|l . being careful to remove 
the unwanted contribution from differentiating the step discontinuity in RU,C{t,t^). 
Using also (|4.9|) and H4.12|l then gives 

2Re = i-dt - 2z{t))mC{t, t^) + d{t - t^)mc{t+,t^) 

- J dt'dt" [{-dt - 2z{t))K{t, t') + 5{t - t')]L{t', t")mC{t", iw) (5.28) 

= - / dt'L{t,t')mC{t' ,t^) + Sit - t^)mC{t+,t^) (5.29) 



5.3. Equilibrium 

Above, we derived general expressions for the energy two-time correlation and 
response, in terms of the known correlation and response functions for the Gaussian 
spins which in turn determine K and L. Before looking at non-equilibrium, we consider 
briefly the equilibrium situation; even here the results for the dynamics are new as far 
as we are aware. 

For the response function one uses that at equilibrium i?r2"C(<) = 
/ (dq) e~^(^°i+")*a;°[T/ (zcq+uj)] for a = 1, 2. Here we have retained a possible nonzero 
equilibrium value Zeq of the Lagrange multiplier, to include the case of equilibrium at 
T ^ Tc- Inserting into H5.27|l and taking LTs gives 



R%\s)^ Hdq) 



where K^q is generalized from H4.22|l to 

k,q{s) =T I (dq) ,..,r„L/., ,,.M (5-31) 



(Zoq + Uj)[s + 2(Zoq + Lo) 

Using the spherical constraint condition T J{dq) (zoq + w)^^ = 1 one then shows, after 
a few lines of algebra, that 

FtJ^(s) = -\ s + 2T- -.J— I = -Lj^J (s) (5.32) 

Remarkably, therefore, the equilibrium energy response function i?^(i) is directly 

proportional to the inverse kernel i}cq{t). Its asymptotics are then given by (|4.23|l . 
The long-time equilibrium susceptibility which encodes the response to a step change 
in the field is, using Kcq{Q) = j [dq) T / [2{zcq + loY] and the generalization of (|4.25|l 

to T 7^ Te, 

X''^ = r dtR%^{t) = 1l(2)(0) = ^ f 1 - TTTTpfT^— ^) (5-33) 
Jo 4 'I 2 V j{dq)[T/{zoq + ujWJ 

It is easily shown that this is consistent with the known result for the temperature 
dependence of the equilibrium energy, E = {H) = J {dq) ^uj[T / {z^q + ui)]: one 
confirms Xe ~ TdE/dT as it should be. The factor T arises because our field h 
is introduced via H ^ H — HH ~ {1 — h)H and so corresponds to a temperature 
change oi T/{1 — h) — T = hT to linear order in h. The inclusion of subleading 
non-Gaussian fluctuations is crucial for achieving this consistency; as discussed at the 
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beginning of Sec. 21 the Gaussian theory does not even give the same answers for the 
fluctuations of H in its two representations as a bond or spin product observable. The 
same phenomenon occurs in a purely static calculation of the energy fluctuations and 
response [T^ . 

The temperature-dependence of the susceptibility H5.33|l deserves some comment. 
As T approaches Tc from above, one has Zoq 0. For d < 4, the denominator of the 
second term in H5.33|l then diverges, and XeI"^ smoothly approaches the value 1/2 
and remains constant for T < Tc- For d > 4, the denominator has a finite limit for 
Zcq 0. This produces the well-known discontinuity in x^/T at T = Tc, since for 
T < Tc the second term in H5.33|l again does not contribute. To see this explicitly, one 
notes that for T < Tc the a; = term in the spherical constraint condition, with its 
weight 1/N, has to be treated separately: 

t( [idq)^— + -^)=l (5.34) 

\J Zcq + UJ NZcq J 

For Zcq (on the scale 0{N'^)) the first integral is j{dq)\/u) = 1/Tc and so 
Zoq = {l/N)TTc/{Tc - T). This then gives for the denominator integral in (jO^ 

which diverges for iV ^ oo at T < Tc as claimed. 

An interesting question is how the discontinuity at T = Tc in c? > 4 of the 
susceptibility x^e^ ^-^^ the integrated response, is related to the temperature variation 
of the time- dependent Ft^[t). In 1)5.3111 one notes that, for T < Tc, Kcq{s) acquires a 
distinct contribution from the w = mode 

kcq{s)= ({dq) ^ +1 ^ (5.36) 

J UJ{S + 2UJ) N Zcq[s + 2Zcq) 

T T —T 

U}{S + 2uj) IcS 

where in the second line we have neglected Zcq = 0{1/N) against s. From (|5.32)l 
one sees that the additional contribution to Kcq{s) produces an extra pole in n^{s) 



UJ 



-2 



which for small (Tc-T)/Tc is located at s = I /to with to = [TTc/4(Tc-T)] / [dq) 
Transforming to the time domain, one finds that this gives an extra contribution of 
~ 1/fo exp(— t/^o) to R^^{t). Crucially, this has a finite integral even in the limit 
T — » Tc, where to diverges, and the appearance of this term causes the discontinuity 
in xl'e Sit T — Tc- Translating these results to the time-dependent susceptibility 

XE\t)= fdt'R'^^it') (5.38) 

one has that, for fixed finite t, x^{t) depends smoothly on temperature around Tc. 
For large t, x^(t) approaches a plateau value equal to the equilibrium susceptibility 
at T = Tc+; the approach to this plateau is as ~ t(4-d)/2^ Yov T < Tc, however, Xb (^) 
eventually increases further on a diverging timescalc t to ^ ^/{Tc — T) to approach 
a higher plateau value given by the susceptibility at T — T^. By FDT (see below) 
the energy correlation function correspondingly shows a power-law decay to a plateau 
for t ^ to, from which it falls to zero only for t > to. 

A final check on our results is that in equilibrium the energy response and 
correlation function should be related by the usual FDT. This is indeed the case. 
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One combines the results ()5.5f) . H5.12|l . H5.14|l . (|5.19|) for the constituent parts of 
the correlation function and decouples all the convolution integrals using temporal 
Fourier transforms. It is easiest to do this starting from expressions where all factors 
of fl have been preserved, e.g. for c'^^ one uses (|5.2|) rather than ()5.5|l . After some 
straightforward but lengthy algebra one then indeed finds the Fourier domain version 
of equilibrium FDT, dtC'j^{t)e-"'* = T[FQ{-iiy) - If^{iv)]/ (iv). 

6. Energy correlation and response: Non-equilibrium, d> A 

We now evaluate the behaviour of the energy correlation and response for the out-of- 
equilibrium dynamics after a quench to criticality. For the correlation function, we 
need the long-time behaviour of 



CnC{t,t^) - / {dq) '^Tliut^f-^ 
J u;^ g 



>3(^w) 



-2ui{t-t„) 



At equilibrium, where Tc — 1 and g{t) ~ const., this gives T^Kcqit — t^). In the 
non-equilibrium case one gets a correction factor which becomes relevant in the aging 
regime, where t — t^j ^ t^j ^ 1: 

CnCit, tw) - TcK,q{t - t^)Tcnc (t/t^) (6.2) 
fty jdwwi^-^y'e-'i^-^^-TUw) 

This is valid for all dimensions d > 2, but in the rest of this section we focus on the 
case d > 4, where k = 0. The regime 2 < d < 4 is more complicated and treated 
separately in the next section. 

Before assembling the four parts of our expression for Csit^t^), it is useful to 
have a guide on what to expect overall. As discussed above, the equilibrium energy 
fluctuations remain finite at criticality. One therefore expects that, in the short-time 
regime t — t^^ CE{t,tvj) will decay as in equilibrium. From l|4.23|l . (|5.32|) and 

FDT it follows that this decay becomes a power law, ^ (t— t^)^'*"'^-'^^ as t—ty, becomes 
large. Aging effects should then appear when t ~ t^ ^ and manifest themselves 
through a scaling function of t/t^. Overall, one expects in the aging regime 

CE{t, iw) = (i - t^)'^^-'^^''^ X [scaling function of t/t^] (6.4) 

and this is indeed what wc will find. Consider now c|j^\ for which wc obtained in 15.5|l 
that 4C^^^ = -{dt + '2z{t))CVtC{t, t„). Now from ijTTI) and the fact that g{t) const, 
for t ^ oo and d > 4 it follows that z{t) — g'{t)/[2g{t)] decreases more quickly than 
1/t. Using also 1)6.2(1 . we see that the z{t) term is negligible for large times and that 

4C^^'' decays as {t — tw)~'^^^ times an aging correction. This is negligible compared 
to (|6.4|l . so we do not need to analyse this contribution further. 
For c]j ' , we have from H5.12|l 

4d^^ = 2z{t)CnC{t,t„) - 2T^CnC{t,t^) + j dt'L^^^t,t')CnC{t',t^) (6.5) 

Because of the factor CQ,C{t' ,t^) it is useful to split the t'-integral into the regimes 
t' > tvj and t' < iw The first regime gives 

Te / dt'L^J'^{t-t')K,q{t' -t„)TL{t/t')Tcnc{t'K) (6.6) 
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As in the analysis of (|4.27|l one can now argue that, for large t — t^, the integral will 
be dominated by the regions t' k, t and t' k, t^. The weights contributed by these 
regions are Kaci{t — tw) /o°° dt'L!^{t') and L'^{t — t„) dt'Keqit'), respectively; the 
integrals are both finite, so that both terms scale as (i — tw)^^"'^''^^- These weights are 
then multiplied by the relevant values of scaling functions in the integrand of H6.6|l . 
The overall result is smaller hy 1 / (t — tw) than the leading term H6.4|l in Ce, and 
can be neglected. The t' < t^ part of the integral from (|6.5|) is 

dt'Lf^it - t')K,^{t^ - t')TL{t/t')Tcnc{t^lt') (6.7) 

The factor Kpq{t^ ~ t') concentrates the weight of the integrand near t' — t^, so that 
for t — tvj and t^ ^ 1 all other factors can be replaced by their values at t' = t^j, giving 
the result TcKcq{0)L^ {t — tyj)!FL{x). This scales as {t — iw)^^"'^-'^^ times a scahng 
function of x and so is also negligible compared to the leading contribution (|6.4|) . 
From H6.2|l the second term on the r.h.s. of (|6.5|l has the same subleading scaling, and 
the first term is even smaller. Overall, Cjj therefore gives a subleading contribution 

to Ce in the aging regime. By very similar arguments one shows that from (|5.14|l 
is also subleading. 

In the aging regime we therefore only need to consider C^j-* from 1)5. 19|) . The 
long-time behaviour of the function CC{t, t^) is easily worked out as 

CC(t,iw) =CCeq(i-tw)^cc(iAw) (6.8) 



(6.9) 



/rTi2 poo 
(dq) ^e-^'^(*-*") = 2T, / dt' K,^{t') ^{t- t^)^^-"^^' 

^cc(x) = ;,^^(.-a)/2,-2(.-,": (6-10) 

where the scaling function Tcci.'x^ ^ x'"^ for x ^ 1. Now consider the f'-integral 
from ()5.19|l . 

I{t, t'J= [ dt' [{2z{t) - 2T,)S{t - t') + L(2) (i, t')]CC{t\ C) (6.11) 



The equilibrium part hf^ {t-t') of L^^) {t^ f) ensures that this has its mass concentrated 
near t' = t. Because we are looking at the aging regime, we have t — t'^{> t — tw) ^ 1, 
and 1)6.9)1 then shows that the function CC{t' , t'^) is slowly varying near t' — t. It can 
thus be replaced by its value there, CC{t, t'^), to give to leading order 

I{tX) = (2z(t)-2Te + L(2)(0))CC(t,4) = -J—CC{tX) (6.12) 

-Keq(U) 

where in the last step we have used z{t) <^ 1. One might suspect that in the t'-integral 
of 1)6. 11)1 . the region around t' k, t'^ also contributes, since this is where CC{t',t'^) is 
largest. However, it can be shown that this region only gives a subleading contribution 
compared to 1)6.12)1 . 

In terms of I{t, t'^) we can write C^"* using ()5.19)l as 

44')-^ / <[(2^(*w)-2re),5(tw-C)+i^'Hiw,C)]^(^,C)(6.13) 
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Again, the integral is dominated by the region t'^ w t^, because of the factor 
Leq {tw — ^w) and we thus get our final expression 

As argued above, in the aging regime the full energy correlator will be identical to 
this. 

We can now turn to the evaluation of the out-of-equilibrium response function. 
We will not write explicitly the second term in H5.29|l . which just removes the unwanted 
S{t — iw) contribution arising from the derivative applied to the first term's step 
discontinuity. Thus, 

2Re = i~dt - 2T^)mC{t,t^) + j dt'L^^\t,t')mC{t',t^) (6.15) 

The long-time scaling of the function RilC is 

mc{t,t^) =mCeq{t~t^)TBncit/t^) (6.16) 

- iw) = Tc J{dq) e-2-(*-*w) ^ Tj^f _ _ _ ^g^^^^ 
_ Jdww('^-^^/^e~^(^-^^-Tc{w) 

The integrand in H6.15|l has its mass concentrated near t' = t, from the factor 
Li^(t — t'), and near t' = t^, from the factor RQ,Ccq{t' — t„). This gives 

2Re = i-dt - 2T,)RnC{t,t^) + L^^^{0)RnC{t,t^) + ii(2)(t - tw)^L(iAw) (6.19) 

using that dt' RQC^qit') = \2T^ dt' J{t') = \K^q{Q) = \ (see after (flTl 'l. 
The last term in 1)6. 19|) scales as (t — tw)^^ '^-'^^ times a function of a; = t/iw, while 
the other terms scale at most as it — iw) ''^^ and so are subdominant. Using also that 

and Leq'' are asymptotically proportional for d > 4 from (|4.24|1 and that = 
then gives the final aging regime expression 

Re = i42)(t-t„)^i(t/t^) = '^L''^\t,t^) (6.20) 

Interestingly, the simple relation H5.32|l between the energy response and i'^^^ therefore 
holds not only at equilibrium but also in the aging regime of the non-equilibrium 
dynamics, and thus overall across the entire long-time regime. 

We can now, finally, find the FDR for the energy correlation and response. In the 
aging regime of interest here, it is useful to rewrite the response function (|6.20() as 

RE{t,t„) = ^ K,q{t~t^)TK{x) = , K{t,t„) (6.21) 
4if2^(0) 4X2^(0) 

using H4.24|l and = J^k- Recalling the definition H4.8|) . and the fact that is 
given by H6.14|l to leading order in the aging limit, we then obtain 

XE{t.,t^)^ „i I . = n^ \r'tr' {b-ii) 
C'E{t,t^) j{dq)C^Cq 

Remarkably, this is in fact identical to the FDR (|3.31|l for large blocks of spin product 
observables. In particular, Xe tends to the limit value X°° = 1/2 for t ^ t^,, and 
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this is identical to the one we found for spin and bond observables: for d > 4, all 
observables we have considered lead to a unique value of = 1/2. 

The FD plot for the energy has a limiting pseudo-equilibrium form; this follows 
from Ht).9|l and Ht).14|l . which show that Csit^t^) has decayed to a small value 
iw"''"'^^ at the point where aging effects begin to appear. More generally, the 
correlation function scales as Ce ^ [t — tw)^'' ''■'^^ for a; w 1, and for a; ^ 1, where 
J'ccix) - a;-2 from (OHIl . as Ce - t^'^-'^^/^x-^ ^ tlt''^^^. 

It should be pointed out that while the FDR for the energy matches that for 
blocks of product observables at long times, the correlation and response functions 
themselves do not. This follows from the nontrivial proportionality factors 1/K^q{0) 
in H6.14|l and H6.21|l . The latter are required to match the limiting behaviour for 
t — ^ of the aging regime results to the asymptotics of the equilibrium results for 
t — <w 3> 1. Indeed, by combining the effective equilibrium behaviour for t — ty, = 0{1) 
with the above aging expressions we can write 

C'Eit, iw) - C^^t - t^)Tcc{x) (6.23) 
RE{t, iw) = R"^{t - t^)TK(x) (6.24) 

and these are now valid throughout the long-time regime, i.e. for large t„ but 
independently of whether t — tv, is large or not. As promised they match with the aging 
expressions for 1 ^ At <C tw For the response this is obvious from H5.32|l . For Ce, 
lIT^ and flFTTHi show that -dAtCj^iAt) = {T^/A)Lf^{M) « [TjAkl^{Q)]Kc^{At) 
for large At; from (|6.9|) this agrees with the corresponding derivative of the 
result in lIHTHl , -dA tCC^^jAt) /[BKl^jQ]] = [2Tjml^{Q)]K^^{At). Note from the 
discussion after (|5.38|l that the function C^'^(At) discontinuously acquires an additive 
constant (non-decaying) part as T crosses Tc from above in d > 4. What we mean 
in H6.23|l is the limiting form for T \ Tc from above, which does not contain this 
plateau. That this is the correct choice can be seen as follows: the non-decaying part 
of C°^{A.t) at T = T~ arises from the fluctuations of the q = Fourier mode, i.e. of 
the magnetization, which are larger by 0{N) than those of the other Fourier modes. 
In the context of our non-equilibrium calculation, all Fourier modes have fluctuations 
of comparable order (in system size), so that this contribution is absent; it would 
appear only for times t^ that scale with system size N. 

Finally, the long-time expressions 1)6. 2316. 24|l also show that an energy FD plot 
would have a pseudo-equilibrium form at long times, because e.g. in C£;(t,<w) the 
equilibrium factor has already decayed to ~ tw"''''^^ of its initial value when aging 
effects appear around Ai ^ 

7. Energy correlation and response: Non-equilibrium, d < A 

In dimension d < 4, the evaluation of the energy correlation function in the aging 
regime is somewhat more complicated. One can nevertheless show that, as before, the 
dominant contribution to is with the other three terms being subleading. 

We thus again need to consider the function I{t' ,t^) defined in (|6.11() . and then work 
out Ce using 16.13|l . One can see clearly that the approach leading above to (|6.12() no 
longer works: in d < 4, Lcq''(0) = 2Tc, so the leading terms in (|6.12|l actually cancel. 
To treat this cancellation accurately, we use (I4.2()|) to write the coefficient 2Tc in the 
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2Tr = -- 
This allows us to express (jfi.llll as 



(7.1) 



2z{t) 
t 



CC{t,t'J 



K,^it) 

+ I dt'Li'^{t-t') 
Jo 

To make progress, we need the behaviour of CC for d < 4. For long times we have 



(7.2) 



:FUt/t')cc{t', C) - ^^cc{t, C) 



2/+ J- \ rTi29{iv 



CCit,t^)= J {dq)C'^{t,t^) = T^ ^^^^ 
The equal-time value thus scales as 



{dq)oj J-Q{u}t^fi)e 



(7.3) 



CC{t,t)^ / da;a;'*/2-3_^^(^t)=t(4-d)/2 / dw w^^-'^^^^T^cM (7.4) 



hence CC{t, t) = -/dt'-^-'^^/^ with some constant ^d- A scaling function is then obtained 
if we normalize CC'{t,t^) by this equal-time value: 

CC{t,t^) _ g{t^) t^^^"^'^ jdww^''-^y^TUyo)e-^^^'-'-y'- 
CC{t,t) git) <(4-'i)/2 Jdww(d-e)/2j:2(^^-) ^ ■ > 

This is a function G{x) of x — t/t^; note that the first two fractions cancel since 
g{t) - t-"" = t('^-4)/2. So far this scaling expression holds for t > t^. To be able to 
use it also for non-ordered times, note that for t <t^, 

=x^''-'^''g{t^lt) (7.6) 



CCit.t^) 
CC{t,t) 
So we have overall 



t^ 



CC{t,t^ 



CC{t,t) 



g{t/t^), g{x) 



t J CC{t^,t^) 

/du;u;(''-6)/2^2(y^)g-2(x-i)» 



for a; > 1 



(7.7) 



for 2; < 1 



One easily works out the asymptotics of Q: Q{x) ^ x ''/^ for x oo, g{x) ^ x" 



^d-2 

for a; ^ 0. For |x - 1| <C 1, on the other hand, 1 - g{x) - |x - l|(4-'i)/2. This 
corresponds to CC{t,t) - CC{t,t„) - |t - t^\^'^-'^^/^ for 1 < |t - tw| < iw (The 
behaviour of this difference for smaller time intervals t — t^ = 0(1) is not captured 
accurately by the scaling form H7.7|) . but integrals over this regime contribute only 
subleading corrections to the results below.) 

We can now insert the scaling form (|7.7|l of CC into (|7.2|l . The non-integral terms 
turn out to be subleading (see below), so 

j(t,C) = J dt' Lf;i{t-t') 

X j^L{t/t')cc{t',t')g{t'/t', 



K,^{t') 



^CC{t,t) / dt'Lf^{t-t') 



t' 



i^cq(i) 



CC{t,t)g{t/t'J 



(7.8) 
(7.9) 
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where we have used that, from and (|OT|) . J^L{t/t')CC{t' ,t')/CC{t,t) = 

(t7t)('*"2)/^(t7t)(**~'')/2 = t'/t. The remaining "subtracted" integral now no longer 
has its mass concentrated near t' = t because the terms in square brackets give a factor 

(2) 

t — t' there. The whole integration range contributes, so that we can replace Lcq {t — t') 
by its asymptotic form L^q^ - t') w A<j(t - t')(''-6)/2 with a d-dependent 

constant A^- Similarly, the ratio Kcqit') / Kcqit) can be replaced by (i'/i)^^"''^/^ to 
leading order. Scaling the integration variable as y = t' /t'^ then gives 



I{tX) = I{x') = ld^dx' 



dy (1 - ylx')(''-''^'^y \g{y) - y-'"^g{x')x""'] (7.10) 



This shows that in the aging regime / depends on x' = t/t'^ only. One can now 
also see that the neglected terms from (17.2(1 are indeed subleading: they scale as 

We can now proceed to simplifying ((6.13() in the aging regime. Using the same 

subtraction method as above, and remembering that Ce ~ C^"^ to leading order, one 
finds by analogy with H7.2|) 



8C£;(t, iw)— / dt'^ L^^Mt-„ — t'^) 



:FL{t^/t'JI{t,t'„) 



(7.11) 



Here subleading terms similar to those in the first line of 1(7. 2|) have already been 
neglected. With the scaled integration variable y = t'^/t^, and using the asymptotic 
forms of Liq{tvj — t'^) and Kcq{t'^) as in H7.1U|I one gets 

SCEit, 4) = A,e-4)/2 dy (1 - y)('^-6)/2 Ud-2)/2 j^^^y^ _ y(2-d)/2^(^)l (7.12) 



This shows that C'e scales as iw ^^-*^^ times a function of x = t/t„. It is difficult to 
work out the whole functional dependence on x. We therefore focus below on the 
asymptotic behaviour for x — > oo, which gives the asymptotic FDR. First, though, we 
check that in the limit a; — > 1, where aging effects are unimportant, our expression 
matches with the equilibrium result C'^{t — t^). From 15.32(1 the latter behaves as 
C^'^(At) = (rc/4) dt' Lf^{t') w [AdTc/2(4- d)]Ai('*-4)/2 ^j. j^^^gg ^^ compare 
with the prediction from ((7.12() . one uses that I{x') ~ ln(a;' — 1) for x' « 1; this 
follows from the behaviour of g{x) for x ~ 1. (Note that the proportionality constant 
in / ~ ln(a;' — 1) is positive^ so that / itself is - in contrast to the case d > 4 - negative.) 
Inserting into (|7.12(l one then finds that for a; « 1 the integral scales as (a; — 1 )('*"'*) Z^. 



Overall, one gets CE{t,t^) ~ t 



(d-4)/2 



(x - l)('^-4)/2 = (t _ i^)(d-4)/2^ consistent with 



the expectation from the equilibrium result. One can work out the prefactor of the 
power law and finds that this, too, agrees as it should. 

Turning now to the behaviour of (|7.12l) for large x, we need the asymptotics 
of I{x'). The leading tail ~ y~'^^'^ of Q{y) is subtracted off in the expression in 
square brackets in 1(7. 10(1 . leaving as the next term y~('='+2)/2. Even together with the 
additional factor y this makes the integral convergent at the upper end for x' — > oo. 
In the limit we thus get I{x') « —ad'yd^dx'~'^ , with 



ad = 



dyy 



9dy 



-d/2 



-Q{y) 



(7.13) 



and gd = lim^ 



> Q{x)x' 



d/2 
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Inserting this inverse-square asymptotic behaviour of I{x') gives for the integral 
in (|7.12|) the scahng a^Pdld^dX"^ with 

(3d = dy (1 - yr-^y^ [i'-'^y^ - i'+^y^) = -tl^J^I^ (7.14) 

It then follows, finally, that CE{t,t^) = {adPdJd^d/^)i'v/^t^'^ for t 3> tw, and 

C'^it, iw) = {adPdld\ld/l6)t^J~^')/H-^ (7.15) 

To get the FDR, we now need the response function Re- This can be evaluated 
very similarly to the case d > 4 and one finds that the last term in ()5.19fl is again the 
dominant one, giving 

i?B(t,iw) = j42q)(t-tw)^L(tAw) = ^L^^Ht^tv.) (7.16) 

(The RflC{t,tyj) terms in (|6.19|l look dangerous: they scale as {t — iw) ''^^ and are 
thus larger than L^^ (t — i„) ^ (t — tw)'''^^''^'^ in d < 3. However, their prefactors 

^ (2) 

— 2Tc + Lcq^(O) cancel. Treating this cancellation more carefully then shows that 
these terms do remain subleading compared to (|7.16|l .l As before, ARe equals i'^' 
in the aging regime, and this then holds across the whole long-time regime since for 
t — tw = C(l) it matches the equihbrium behaviour ARE{t,t^) — Lcq {t — t^j). For 
t ^ tw, on the other hand, the response (|7.16() becomes 



/ . \ (2-d)/2 , 

RE{t, t^) = \\dt^'-'^'' j = 



2 



(7.17) 



Comparing with H7.15|l then finally gives for the asymptotic FDR for the energy in 
d<4, 

AT 

= . \ , (7.18) 
OidPdJd'^^dd 

After evaluating the various numerical factors (see [Appendix Aj l this can be written 
as 

dad T {^) 

where 



= IT- ' \.L\ ' (7-19) 



OLd 



■ . . / dx fdyy^'^-'^'^ ^'^-f""' (l - , - .-('^+^)/^) (7.20) 

sm[7r(4 - d)/2J Ji Jq 1 + x-y \ I 



Near c? = 4, one can expand in e = 4 — d. It can be shown by explicit calculation that 
the integral in 17.20|l is exactly zero for d = 4, giving ctd = 2/e + 0{e) and so 

^ ;r(4-6)r(e/2) = i _ 1 + o{e') (7.21) 

^ (4-e)r(4-e/2) ^'^23^'^ ^ ' 

Notably, this is different from the FDR = 1 - 2/d = 1/2 - e/8 + 0{e^) for aU the 
other, finite-range, observables that we considered previously in d < 4. It is, however, 
consistent with RG calculations to 0{e) for the 0{n) model in the limit n ^ oo, 
for an analogous choice of observable [^. Remarkably, therefore, the non-Gaussian 
fluctuations induced by the weak infinite-range interaction in the spherical model seem 
to mimick precisely the effects that are seen in more realistic models such as the 0(n), 
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Figure 3. Asymptotic FDR X°° vs dimension d for finite-range observables 
(equation lij.3l . dashed) and energy (equation 17.201 . solid). Dotted lines indicate 
the corresponding linear expansions 1/2 — e/8 and 1/2 — e/3ine = 4 — d. 

even though in the latter aU interactions are short-ranged and there is no difference 
between the behaviour of block observables and global ones. 

We show in Fig. |3| the dependence of the asymptotic energy-FDR on 
dimension d and compare with the result for finite-range observables. They agree 
in d > 4, but the difference between them grows as d decreases towards d = 2, with 
the energy FDR always having the lower value. In the limit d ^ 2, both FDRs 
converge to zero, but while the finite-range FDR X°° — e'/2 -I- 0{e'^) does so linearly 
in e' = d — 2, the energy FDR vanishes quadratically as — + ©(e''^), due to 
the divergence of ad = + ©(l/e'). 

As in the case d > 4, an energy FD plot would have a pseudo-equilibrium form 
which hides all non-equilibrium effects at long times. Indeed, one could write (|7.12|l 
in the form CE{t,t^) = C'^{t — tw).?^Cf; (a^), and the decay of the equihbrium 
factor C'^ squeezes all details about the aging factor Tce{^) into a vanishingly 
small region of the FD plot for long times. While we have not calculated !Fce 
explicitly, the discussion after H7.12|l shows that TceO-) = 1 as it should be. For 
large x, on the other hand, we saw in (|7.15|l that CEit,tw) ~ ti/^t~'^, implying 
^Ceix) - tt!'^t-'^{t - iw)^""''^/^ ~ x-'^^^. This matches continuously at d = 4 with 
the prediction l().23l) for d > 4, where the aging correction decays as Tcc{x) ^ x^"^ . 

8. Magnetized initial states 

In this final section we consider the dynamics for initial configurations with nonzero 
magnetization, focussing as before on the non-equilibrium dynamics that results when 
the system is subsequently brought to the critical temperature T^. Physically, such a 
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situation could arise in an "up-quench" , where the system is equihbrated at T < 
and temperature is then increased to T = Tc. As explained in the introduction, our 
interest in this scenario arises from recent results [30] which show that such initial 
conditions produce FDRs that differ nontrivially from those for zero magnetization. 
The analysis of [201 was limited to high dimensions d or infinite-range interactions, 
however; the calculation below will allow us to see explicitly how the results change in 
finite dimension. In particular, we will obtain exact FDRs for magnetized coarsening 
below the upper critical dimension, i.e. d < 4 in the spherical model. 

We will continue to use the notation Cq(i,tw) = i^/^) {Sq{t)S^{t)) for 
Fourier mode correlations. For q = this is now a full, unsubtracted correlator, 
with Co{t,t) = Nm'^it) + 0(1) and m{t) = {l/N){So{t)) the time-dependent 
magnetization. The difference Co(t, iw) — Co{t,t^) — Nm(t)m{t^) is the connected 
correlator, which has values of 0(1) and is the relevant quantity for analysing the FD 
behaviour. For q 7^ 0, connected and full correlators coincide: 



C'q(t,iw) = Cq(t,iw) - [l/N) {S^{t)) {S*Jt^)) = C^{t,t^) 



(8.1) 



We now need to check how the analysis in the previous sections is modified by 
the presence of a nonzero magnetization. The Fourier space equation of motion H2.5|l 
remains valid, and so do the resulting expressions for the response function 
-Rq(i,tw) (ESI) and the fuU correlator Cq(i,iw) H2.10I2.13|I . The expression (|2.16|l 
for the Laplace transform of g{t) that results from the spherical constraint also still 
holds, but the solution is now different. In the q-intcgral, the q = contribution 
(l/A^)Co(0, 0)/s = m^(0)/s has to be treated separately. In fact, one sees that for 
s ^ this term always dominates the rest of the integral, which diverges less strongly. 
At criticality, where T = Tc = [/ {dq) l/uj]-^ = [2/(0)]"\ one thus has for s -> 



2T,g{s) 



{dq) 



1 

2cj 



1 



2uj 



which using H4.22(l can be rearranged into 



m; 



'(0) 



[dq) 



LU{2U! + s) 



cq 



■2) 



(8.3) 



For d > 4, i^eq(O) is finite so this scales as s^^; for d < 4, on the other hand, Kcq{s) 
diverges as s(''~4)/2 so that g{s) ~ s~'^/^. Translating back to the time domain, g{t) 
behaves for large t as 



g{t) - m^{0) i". 



(8.4) 



1 (rf > 4) 

{d-2)/2 (d<4) 

Note that this asymptotic behaviour is independent of any details of the initial 
condition except for the presence of a nonzero m{0); it depends on the actual value 
of m(0) only through the prefactor m^{Q). For the time-dependence of m{t), one gets 
by taking an average of (|2.8|l 

m(0) 



m(t) = i?o(i,0)m(0) = 



(8.5) 



Because of the proportionality of g{t) to m^(0) for large t, the asymptotic decay of 
m{t) is independent of the initial conditions, in terms of both the decay exponent and 
the prefactor. 
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8.1. Finite-range spin observables 

We first analyse the correlation and response functions for observables that relate to 
a number of spins that is much smaller than N. As for the unmagnetized case, the 
fluctuations of the Lagrange multiplier z can then be neglected. To understand the 
magnetized case, it is useful to shift the spin variables by m{t). We will see that the 
equations of motion then acquire the same form as before, so that we can directly 
transfer the main results from the unmagnetized case. Explicitly, we consider the 
following decomposition of the spin variables 

S,=m + [/, (8.6) 

where Ui is a zero-mean variable. The equation of motion (|2.4|l for Si then gives 

dtm + dtUi = {m + Uj)+ - z{m-\-Ui) (8.7) 

From (IH31) and the definition dt\inn{t) = -{l/2)dth-ig(t) = ~z{t), so dtm = 

—zm. Also ^ij — Oi giving 

dtU, = -n^,U, + - zU, (8.8) 

This is the same as the equation for Si in the unmagnetized case, and so one can deduce 
directly the solutions for the correlation functions of the Ui; these are the connected 
correlations C. The initial values again become unimportant for long times, allowing 
us to work out the scaling of the C, and then together with the response R also the 
FDR X. It is clear from the description in terms of the subtracted spins Ui that there 
is nothing special about the case q = 0, and all results will have a smooth limit as 
q ^ 0. Because we are neglecting the fluctuations of the Lagrange multiplier z, this 
limit again has to be understood as that of a block magnetization calculated over a 
region much larger than the correlation length (in the time regime being considered) 
but much smaller than the linear system size, so that 1/L. 

Applying fOU) . we can now write down directly the connected correlation 
function at equal times as 

C'q(iw,tw) = 6-^(0, 0)— — + 2Te—— / dt'e^^''g{t') (8.9) 

9[tw) g{t„) Jo 



At long times, the first term is subleading due to the scaling H8.4|l . and one has the 
behaviour 

C'q(iw,iw) = — ^c(c^tw), Tc{w)^2w[ dyy"e-2-(i~?/) (8.10) 
^ Jo 

This result is of course the same as ()3.4f) . except for the replacement of by a which 

reflects the difference in the asymptotic behaviour of g{t). 

The two-time connected correlations are Cq{t,t^) = Rq{t,t^)Cq{t„,t„) with 

-Rq(i, iw) given by H2.6|l as before. As a consequence, the expression H3.1|) for the FDR 

Xq also remains valid, and one finds the scaling form Xq{t, i„) — Txi^t^) with 

Tx\w) = 2-{2w + a)[ y^e-^""!!-!/) (8.11) 

^0 

which is directly analogous to (|3.4|l . In the limit w = ujt„ oo, — 1, which means 
that all modes with q ^ equilibrate once t„ ^ 1/lo. In the opposite limit w — > 0, 
corresponding to t„ ^ l/w, 

Xq.^^ = | 2/3.,^^ (^>f) ^8.12) 



d/{d + 2) (d<4) 
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For q ^ this resuh appUes independently of the value of t^ as long as 3> 1, so 
that the FDR for the block magnetization will be a straight line with slope (|8.12|) . This 
is as for the unmagnetized case, but the actual value of the FDR is now different. It is 
also different from the value X°° — 4/5 predicted for Ising models in the limit of large 
d [201; we will see below that the latter value is obtained for the global magnetization, 
which is affected by local spin fluctuations of 0{N~^/'^). 

For later reference we write down the long-time forms of the correlation and 
response functions for q ^ 0. By setting u; = in 1)8.9(1 and taking the long-time 
limit where the first term becomes negligible, we find for the connected equal-time 
correlator 

C'o(tw,tw) = ^ (8.13) 
1 -I- a 

The response function is, from 1(2.6(1 and 18.5(1 . 



where the last equality holds for long times. The two-time correlator is therefore 

2Tt ft ^ 



Co(t,iw) = ^(^f j (8.15) 

From these results one of course retrieves the long-time FDR Xo(t,tw) = 
Tci?o(t,iw)/C'o(t,tw) = {a + l)/{a + 2), obtained in via the limit ujt^ -> 0. 

As explained above, these results apply in the regime 1/i <C g ^ 1. For q = 
itself, they capture only the Gaussian part of the spin-fluctuations, and non-Gaussian 
corrections become relevant as discussed in the next section. 



8.2. General expressions for magnetization correlation and response 

We now turn to the FD behaviour of the global magnetization, corresponding to q = 
rather than 1/L. All N spins are now involved and one needs to account for the 
fluctuating contribution of the Lagrange multiplier, which we write as z + N^^^^Az as 
before. To understand why this is necessary in the magnetized case, but was not in the 
unmagnetized scenario, consider the equation of motion 12.511 for the zero-wavevector 
Fourier component of the spins, 

dtSo = -{z + N-^/^Az)So+^o (8.16) 

In the unmagnetized CRS6, 5*0 is cL zero-mecLii quantity of ©(TV^/^). The Az-term then 
contributes only subleading 0(1) fluctuations. For nonzero magnetization, on the 
other hand, the mean of So is Nm, with fluctuations around this of 0{N^/^). The 
coupling of Az with m then gives an 0{N^^^) contribution to dtSo, which is no longer 
negligible. 

To find the resulting non-Gaussian fluctuations in Sq explicitly, we make the 
decomposition Si ~ Si + N^^^'^r^ as before. The discussion in Sec. 01 then goes 
through, and we retrieve 1(4.11(1 for the C'(iV~^/^)-corrections to the spins. For 
the zero Fourier mode, in particular, we have 

ro{t) = -\j dt'df'Roit, t')so{t')L{t\ t")A{t") (8.17) 

To simplify the calculation of connected correlations, we now decompose the Gaussian 
part of the spins into Si = m + Ui, so that the Ui are zero-mean Gaussian variables. 
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This corresponds to a decomposition [/,; — Ui + y/ Nri of the fluctuating parts of the 
spins into leading Gaussian terms and small non-Gaussian corrections, in analogy 
to the representation 5^ = 5^-1- VNvi in the unmagnetized case. The Ui obey the 
equation of motion (|8.8|l , and their correlation and response functions are the Cq and 
/iq calculated previously. 

We will write the connected correlation function for the global magnetization 
which includes non-Gaussian corrections as Cin(t,tw) = (1/-^) {Uo{t)Uo{tvi)) ■ Making 
the decomposition into Gaussian and non-Gaussian parts, this reads 

a.(i,tw) = ^{[uo{t) + N-^/^ro{t)][uo{t^) + N-^/^ro{t^)]) (8.18) 



with, from H8.17|l . 



2 

N 

Here we have defined 



ro{t) ^ \ i dt' dt" Ro{t, t')[Nm{t') + uo{t')]L{t' , t")A{t") (8.19) 



dt' M{t,t')A{t') (8.20) 



M{t,t^)^ j dt' Ro[t,t')m{t')L{t' ,t^) ^ m{t) J dt'L{t\t^) (8.21) 

where the second form follows from (|8.14|) . M is, like R and L, causal and so vanishes 
for t < t^. In H8.20|l we have also discarded the Gaussian fluctuation term uq, which 
is of 0{N^^'^) and so negligible against the term Nm. This is in line with the intuition 
discussed earlier that non-Gaussian fluctuations arise only from the coupling of Az to 
m. Note also that ro is 0{N), so that in (|8.18() the non-Gaussian correction N~^^'^ro 
is of the same order as the Gaussian fluctuation uo, again as expected. 

Substituting H8.20(l into H8.18(l . we see that we need the two-time correlations of 
Uo and A. In the presence of a nonzero m the latter becomes 



/N ^ VN ■ 

The required correlations are therefore {uo{t)uo{t^)) = NCQ{t,ty^) and 

(uo(t)A(i')) = J2^u.,{t)[m^{t') - 1 + 2u^{t')m{t') + u){t')]) (8.23) 
2 

= —m[t')^{u,{t)uj[t')) = 2m{t')y/NCo{t,t') (8.24) 

For the autocorrelation of A, we can exploit the fact that (A) = to write H8.22|l as 
A = ^.(2M,m + ul -{uf)). This gives 

(A(i')A(C)) = 1 ^([2«,(t')m(t') + uj{t') - (u?(t'))][- • • - C • • •]) (8-25) 

ik 

= ±m{t')mit'J Y.{u,it')u,it'J) + A ^(^,(i')«,(C))' (8.26) 

ij ij 

= Amit')mitl)Co{t',t'J + 2 f {dq) Clit' ,t'^) (8.27) 
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where we have used Wick's theorem to simphfy the fourth-order average {{uf — 
~ (^p)) — i^i"^^) ^ — '2{uiUj)'^. Abbreviating the g-integral as 

CC{t' ,t'^), the full connected correlation function (|8.18l) can thus be written as 

an(i,iw) = Co(i,tw)--^ / dt' M{t,t'){uo{t^)A{t')) 
2v N J 

1 



2VN 



dt' Mit^,t'){uoit)A{t')) 



+ i / di'<M(t,t')M(tw,C)(A(C)A(i')) (8.28) 



4 

^C^\t,t^) + C^\t,t„) (8.29) 

where 

CW(t,iw) =Co(i,tw)- / dt' [Mit,t')Co{t^,t') + Mit^,t')Coit,t')]mit') 



dt' dt'^ M{t, t')M{t^, C)m(t')m(C)Co(t', t'J (8.30) 

dt' dt'^ [5{t - t') - M{t, t')m{t')] 
X [Sit„ - t'J - Mit^,t'Jm,{t'J]Co{t',t'J (8.31) 

and 

Ci\t,tJ^]^j dt'dt'^M{t,t')M{t^,t'JCC{t',t'J (8.32) 

Next we derive an expression for the corresponding magnetization response 
function. To this purpose we expand the spins for small fields as 

S, = s, + hr, (8.33) 

where Si are the unperturbed spins and we neglect the 0(iV~^/^) non-Gaussian 
corrections as irrelevant, as in the unmagnetized case. The Lagrange multiplier is 
similarly written as z -I- HAz. By collecting the 0{h) terms from the equation of 
motion for the Si, we then find by analogy with (|5.21|) that the obey 

dtri — —ilijrj — zri — Az Si -\- 5{t — tj (8.34) 

Here the last term represents a field impulse at time t^, uniform over all sites as is 
appropriate for the field conjugate to the global magnetization. Since ri{t) = before 
the field is applied, i.e. for t < t^/, this impulse perturbation gives ri{tj = 1. Starting 
from this value we can then integrate (|8.34() forward in time to get 

r,(t) = Vi?,j(t,tw)- / dt' R,j{t,t')Az{t')sj{t') (8.35) 

The condition we need to impose in order to get Az is that the spherical constraint 
(1/iV) + ^^i)^ = 1 needs to be satisfied to linear order in h, giving the condition 

il/N)Y,^{r^s^) = 0. Inserting into this yields 

Ro{t,tJm{t) ^ dt' K{t,t')Az{t') (8.36) 

where we have used the definition H4.8|) of K{t, t'). Applying the inverse kernel L gives 

Az{t)^ [ dt' L{t,t')m{t')RQ{t',tJ (8.37) 
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Note that this resuh vanishes when m = 0, consistent with the fact that we did not 
need to consider perturbations of z in our calculation of the magnetization response in 
the unmagnetized case. We can now write down the magnetization response function, 
which we denote by Rmit, tw)- It is given by i?m = (1/-^) 'l^ii^i)' inserting the result 
for Az into (|8.35(l . we get explicitly 

Emit, tw) = Ro{t, t„) - J dt" dt' Roit, t")m{t")L{t", t')m{t')Ro{t' , t^) (8.38) 

= J dt'[6{t-t')- M{t,t')m{t')]Ro{t\t^) (8.39) 

This completes the derivation of the general expressions for the magnetization 
correlation and response. To make progress, we need to find the kernel M{t,t'). 
This requires L{t,t'), which is the inverse of K{t,t') = /(dg) i?q(t, t')Cq(t, i'). As is 
clear from the discussion in Sec. 01 the correlator occurring here is the unsubtracted 
one. Because Co{t,t') = Nm{t)m{t') is 0{N), the q = term needs to be 
treated separately in spite of its vanishing weight It makes a contribution 

{l/N)Ro(t,t')Nm{t)m{t') = ■m'^{t)9{t - t'), where we have simplified using (|8.14() . 
We can thus write 

K{t, t') = k{t, t') + n?{t)e{t - t') (8.40) 

with the q 7^ contribution 

k{t, t') = j (dq) i?q(t, t')Cqit, t') (8.41) 

We have switched to the connected correlator here; this makes no difference for q ^ 0, 
but allows us to include q = in the integral because Co — 0(1). To say more, we 
will need to distinguish between dimensions d> A and d < A. 



8.3. Magnetization correlation and response: Non- equilibrium, d > 4 

The scaling of the connected part K oi K can be analysed exactly as in the case of 
zero magnetization: it consists of the same equilibrium time dependence modulated 
by an aging function, K{t,t') = Keq{t — t')T j^{t/t'). The aging part can be worked 
out exactly as in H4.15(l with the only difference arising from the changed asymptotic 
behaviour of g{t) ^ t" rather than i"". For JF^ we can therefore use directly (|4.1ti|) . 
with K replaced by —a: 

2 

In d > 4, where a = 1, the integral can be computed explicitly to give 

d~2 /x-l\ 2 /"x-l 
d-A \ X ) ^ d- 4 



jr^(a;) ^ __(a;- l)(''-2)/2^-a ^ dy y^{x ~ y)-'^^^ (8.42) 







•^k(^) -1-^(^)+T^(^) (8-43) 



We will see below that the precise behaviour of this function docs not affect the 
results. Briefly though, for 2: — 1 ^ 1 the second term on the r.h.s. is leading so that 
J^f^ decreases linearly with a; — 1, while for large x one finds by expanding in l/x that 
^^«[(d-2)/4] 



-2 



To find the inverse kernel L, consider how K{t,t') varies with t. The first part 
in H8.4Q|I starts off close to unity and decays on 0(1) timescales t — t' as Kcq{t — t'), 
with a modulation by the aging factor !Ff^{t/t') once t — t' becomes comparable to 
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t' . The second part, on the other hand, is smaU initiaUy but only decays on aging 
timescales. Comparing K^qit — t') ^ {t — t')(2-<i)/2 to m^{t) ^ 1/t, this second term 
therefore eventuaUy becomes dominant, for t — t' ^ ^2/(^-2)^ This discussion suggests 
that also the inverse kernel L should be composed of two parts with distinct long-time 
behaviour. We therefore write 

L^L + Lo (8.44) 
where L is the inverse of K and Lq arises from the zero-wavevector part of K . The 
continuous part of L is then similarly decomposed as L^^^ = L^^^ + ^q'^''. 

We proceed by writing the defining equations for L^^^ and L'^^\ The full inverse L 
is defined by (|4.9|l and as before has singular parts which are related to the behaviour 
of K{t, t') for t — > t' . One can show directly from the definition of K, and exactly as 
in the unmagnetized case, that 

K{t'+, t') = 1, dt.K{t, tX^v^ = 2Te (8.45) 

The decomposition (|4.12(l of the inverse kernel L therefore also remains valid, and 
from H4.9(l and H8.4()(l we get the following equation for its continuous part L^^' 

dt'[k{t,t')+m^{t)]L'^^\t' ,t^) = 2Tei^(t,tw) +2Tem2(t) - dt,K{t,t^) (8.46) 

This is the analogue of the relation (|4.19|) for the case m — 0. We can argue similarly 
for L, which is defined by 

dt' K{t, t')L{t', tw) = 5{t - t„) (8.47) 

From l|On|l and 

K(t'+,t') ^l-m^U), dt'K(t,t') ^2Tc (8.48) 
and this initial behaviour implies that L can be decomposed as 

l-m2(iw) [l-m2(t 
Inserting into the definition (|8.47|) gives for the continuous part Z^^^ 

2T 1 

dt' K{t,t')L^^\t',t^) = _ / ..^ K{t,t^) - -—^—dt„k(t,t^) (8.50) 

Now for long times, we can approximate l—rrP'^t^) k, 1. Then H8.5()|l becomes identical 
to the relation H4.19|l which determined i'^^-' in the unmagnetized case. Since K has 
the same scaling form as K in (|4.19(l . except for the replacement of Tk by JF^, the 
solution for Z^^^ can be found in exactly the same way. In particular, the scaling 
functions describing the aging corrections in K and L'^^^ are again identical, and we 
can write directly 

n^\t,t') = Lf^{t-t')T^{t/t') (8.51) 

as the long-time form of L^'^\ Here icq-* is the same function as in the unmagnetized 
case, with Laplace transform (|4.21l) . 

It now remains to find Lq^-* . Subtracting H8.50|l from H8.46|l gives 

dt'k{t,t')L^^\t',t^)+m\t) f dt'L^^\t',t^) = 2y^ -2"|^(^w) + m^(^w) ^(^^^^^ 
+ m\t) [2Te - /* dt'L(^\t',tA + -^^^%^9t„K(t,tw) (8.52) 



j^C'.'w) = . + u"L w ''<'' - '») - i"'(i'.<w) (8.49) 



/' 
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To make progress we assume that, by analogy with the structure (|8.40ll of K, Lq^' 
varies only on aging timescales; we will find this confirmed a posteriori. We can 
then concentrate on the aging regime t — iw ^ 3> 1. In this regime, the integral 
/j* dt' W^\t',t^) = /j* dt' Li^^it' - t„)T^{t'/t^) on the r.h.s. of l|H3a becomes to 

leading order J^"^ dt' Lf^{t' - t^) = (0); the aging correction J^f^ is unimportant 
because the integral converges for t' — t„ ^ 0{1), and the upper integration limit can 

" ('21 

be set to infinity for the same reason. From (jOT|l . L^q^(O) = 2Tc - and 
so the square bracket on the r.h.s. of H8.52I) becomes simply the constant 1/Keq(fi)- 
In the first and third term on the r.h.s., on the other hand, K and dt^K scale as 
{t — ty,)~^'^~^'>^^ and {t — tw) '^^^i respectively, and are negligible compared to the 
second term in the aging regime. Disregarding these subleading terms, equation (|8.52|l 
is transformed to 

dt'Kit,t')L^^\t',t^) + m^it) f dt'L'i\t',t„)^^I^ (8.53) 

(2) 

In the first integral, if as assumed Lq varies only on aging timescales, the integral 
is dominated by the region t' t because of the factor Koq{t — t') in K{t,t'). This 
factor again makes the integral convergent within a region t — t' — 0{1), and so we 
can approximate it by Kcq{0)LQ '{t, iw)- This gives 

dt'L^^\t',t^) = ^ - ^^oHt,t^) (8.54) 

and specifically in the limit t/tw 1 

4'^(*w,tw) - K-^\0)m\t^) (8.55) 

To find Lq'^'' for t/tw > 1, we use that m?{t) = jia/t for large t and in d > 4, 
see (|8. 418. 511 . with some dimension-dependent coefficient /i^. Taking a derivative 
of H8.54fl with respect to t then gives 

L^o\t,t^) [l + Md'-ft:cq(0)] = -tnd-^k,^{Q)dtL'i\t,t^) (8.56) 
This implies dilnh'^Q^) / d{\nt) = -[1 + Aid/i^cq(0)] and so together with 

L^^\t,t^)=k-^{0)1^(±] (8.57) 

This can be simplified because in fact fjid = ^oq(O). To see this, note from (|8.3|l that 
g{s)/m?{Q) = l/[i^oq(0)s^] for small s; here we have used that i^oq(O) is finite for 
d> A. Transforming back to the time domain gives g{t)/m'^{0) — m^'^{t) = i/Xcq(0) 
for large i, i.e. m'^{t) = Xcq(0)/t. Thus (|8.57|) simphfies to 

4'^(i,M = ^ (8.58) 

(2) 

This result is consistent with our assumption that Lq only varies on aging timescales. 
Overall, we have thus found for L^'^'>{t^t^) the following long-time form 

(i, t^) = [t - iw)^K (V^w) + (8.59) 
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Of course, for time differences t — — 0(1), Lq {t,t^) will deviate from the 
form H8.58|l derived for the aging regime. However, one can verify by expanding 

both sides of (|8.52|) to 0(i — f^) that, even for t = f^, l'^q^ remains of order 1/iw, so 
that these small deviations are always negligible compared to the first term in (|8.59(l . 

We can now proceed to find M as defined in H8.21|l . and from there finally the 
explicit forms of the magnetization correlation and response functions. Using the 
general structure (|4.12|1 of L we have 

m-\t)M{t,t^)=6{t-t^) + 2T^~ f dt' L^'^\t',t„) (8.60) 



m 



The integral can be separated into the contributions from the two parts of i*^^-* as 
given in (|8.59|) . The first part yields an integral that converges ioi t' — = 0{1); for 
t — 3> 1, it therefore gives Lcq{0) = 2Tc — l/7^eq(0) = 2Tc — l//id to leading order. 
The second part, on the other hand, yields expHcitly dt' {tv; / fJ-dt'^) = fJ,'^^ {l~tw/t), 
so that 

^\t)M{t,U)^5{t-t„) + — -—(l-^-^] =S{t-t„) + —^-^ (8.61) 

fM fid \ t J t 

This result applies for t — t^, ^ 1. For t ~ t^, — C>{\) it is not accurate; e.g. at 
t = tv, the continuous part of ■m^^{t)M{t,tv,) is, from H8.6()|l . 2Tc rather than l/fid- 
However, this deviation over an 0(1) time-range only gives subleading corrections in 
the integrals over M that we need below, as indeed does the 6{t — t^)-term. This can 
be seen in (|8.31l) and H8.39|l . where only the combination S{t — ty^) — M{t,t^)m{tyj) 
occurs; the latter can be written as 

1/2 1/2 

5{t - iw) - M{t,t^)m{K) = 5{t - iw)[l - m(i)m(tw)] ~ ^ ^ Kt - ^w) - ^(8-62) 

The second form holds for long times, where the m{t)m(t^) term that originated from 
the (5-term in H8.61|l is negligible. 

We can now work out the expression H8.29|l for the full connected correlation 



function. One can show that the contribution Cm involving CC is negligible in 
the long-time limit. In the expression H8.3()|l for the remainder Cm"*, let us call 
the second and third term I2 and I3. We need Co, which from (|8.15|l reads 
Co(t',t'^) ~ Tc^w^^t'"^/^ for t' > t'^; because Co is symmetric in time this then 
implies Co{t' ,t'^) — Tct'^^'^t'^, for t' < t'^. Paying due attention to this temporal 
ordering of the arguments of Co and using H8.62|l . one finds 



j.n/2 J./3/2 ft ,11/2 j.3/2 

dt'L^' + / dt' " 



^3/2 1/2 / ^3/2 ^/l/2 



^/l/2 ^/3/2 

^c/ dt'-^-^ (8.63) 



fA. 2. 



Similarly, the double integral in (|8.3Q(I can be evaluated as 

pt i-t, ./1/2/1/2 
/3(t,tw)= / dt' d4_^Co(C,0 (8.65) 



Jo " tv^ tT 

ft -i/l/2 j-'l/2 ,'3/2 ^t' ,/l/2 /1/2 ,'3/2 

T.I dt' i <^^^ + 2Tey^ dt'J^dt'J-^'-^'-^ (8.66) 
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f(T)" r^-'i^w 

3/2 

(2t~tw) (8.68) 



Our final long-time result for the connected magnetization correlator including non- 
Gaussian corrections is then 

3 /2 

an(i,tw) - +/2(t,iw) +/3(i,tw) (8.69) 

f \ •^/^ / 4 2 1 1 \ Tf /t ^^^^ 

Lvir \ / t:. -L -L \ J-pL-Mir / I'll 



For the conjugate magnetization response function we get from (j8.39|l and (|8.61(l . 
together with Ro{t,t^) = (tw/t)^^^, 

Rrnit,t^)^ Roit,t^)^ J dt'_i?o(<',iw) (8.71) 

xl/2 1/2 / ^3/2 



^(^-M-^fj (8.72) 
The FDR follows finally as 

Interestingly, this agrees exactly with the result for the Ising ferromagnet in the limit 
of large dimensionality d |^. As in the unmagnetized case, we see therefore that it 
is the global observables in the spherical model, which are strongly affected by non- 
Gaussian fluctuations, that behave like their analogues in short-range models. In 
fact, even the expressions for the correlation and response functions we find here are 
identical to those in the large-d Ising case, implying "universality" at a more detailed 
level than one might have expected. 

It is worth noting that the effect of the non-Gaussian corrections is very large: 
compared to the Gaussian result Co{t,t^) ^ t„{t^/ty/^, the corrections increase the 
decay exponent to Cm{t,tvj) ^ tw(twA)^^^, so that Cm/Co ~ ty,/t <C 1 for t » ty,: 
there is an almost perfect cancellation of Gaussian terms and non-Gaussian corrections 
for well-separated times. Similar comments apply to the response. The overall effect 
of the non-Gaussian corrections on the FD relation is to leave this as a straight line 
(since X is constant), but to increase the slope from 2/3 to 4/5. 



8.4- Magnetization correlation and response: Non-equilibrium, d < A 

We now consider systems below the upper critical dimension, o? < 4; here there are 
no predictions yet from other models for the non-equilibrium FD behaviour following 
a quench of a magnetized initial state to Tc (but see Sec. 01. As in the case of the 
energy correlations for unmagnetized initial states, leading order cancellation effects 
have to be taken care of for these low values of d. 

We again need to know the scaling of K to determine L^^^; from this we then 
get M and finally the correlation and response functions. The long-time scaling of 
the connected part of K is K{t, t') = Kcq{t — t')Tf^(t/t') as before, with Tp^ given by 
H8.42|l but now a — {d — 2)/2. The contribution to K from q = 0, given by second 
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term in H8.40|l . is negligible relative to k{t, t') for t~t' = 0{1). However, for t-i' > 1, 
it becomes comparable and has the same overall time scaling as K{t,t') in the aging 
regime. To see this explicitly, recall from (|8.4I8.5|) that the square of the magnetization 
decays asymptotically as m?{t) = ^J.dt^°' — fJ,dt^'^~'^^^'^ with some constant /i^- 
Similarly, the equilibrium part of K behaves as Kcq{t — t') — kd{t — i')(2-'')/2 fQj- 
t - > 1 (see after H4.13ll V This gives 

m\t) ^ ^,,[t-t'r-^y^ ^ ^ ( tit'-i \ 

Xeq(t - t') kdtid-^)/^ kd \ t/t' J ^ ■ ' 

in the aging regime where both t' and t—t' are large. (For t—t' ~ 0{1) this expression is 
not accurate but this is irrelevant because there the term (t) is subleading compared 
to K{t,t') anyway.) We thus have the overall scaling of K 

K{t,t')^K,q{t~t')TK{t/t'), TK{x)=T^{x) + t^(^\ (8.75) 

To simplify the scaling function, we integrate the expression (|8.42|l by parts and rescale 
y — > xy, bearing in mind that a = {d — 2)/2: 

(d-2)/2 r ^ 

.(^_l)(2-<i)/2 



d~2 / x-1 



X 

l/x 



d-2 







(8.76) 



For X 1 the integral becomes a Beta-function which evaluates to T{{d— 2)/2)r((4- 
d)/2), and extracting this term gives 

(d-2)/2 

^^(a;)=x(2-'')/2-r(|)r(4^^ 



2 

d-2 /x-1 



1 



X 

r-^\ / dyy^'-^^l\l-y)^^''^l^ (8.77) 

The second term has the same dependence on x as the additional contribution from 
zero wavevector in (|8.75ll . In fact, it turns out that these two terms cancel exactly: 
from our definitions of k^ and [id we have Kcq{t) = kdt'^'^~'^^^'^ for large while 
m~^{t) = g{t)/m'^{0) — [i^^t'^''-^^^/^ . Laplace transforming gives i^eq(s) = fcdr((4 — 
d)/2)s(''~'')/2 and g{s)/m'^{0) ^ ^^^r(d/2)s-'^/2 to leading order for smaU s. But 
then shows that ^;^^r(d/2) = [kdT{{A-d)/2)]~'^, or Hd/kd = r(d/2)r((4-d)/2), 
proving the cancellation anticipated above. Overall, we thus have for the scaling 
function of K 

:Fk{x) = + 'LA j j^^ dyy^'-^y^l - V^-'^^' (8-78) 

Expanding for a; w 1 , one sees that the leading order variation is linear in x — 1 : 

w) 1 + ^(1- - 1) + ^ " 

Note that the prefactor is positive, so that !Fk{x) increases with x in the current 
scenario. This trend persists for all x, not just x « 1, and the scaling function 
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monotonically approaches a limit value for x — > oo. The latter follows from H8.75|l as 
l-^d/kd ~ r((i/2)r((4 — d)/2), since the connected contribution T^{x) decays to zero 
for X oo. 

With the single overall scaling (|8.75|) of K we no longer need to decompose L*^^^ 
into L^'^^ and L^q^ as we did in d > 4; instead the long-time behaviour of L^^^ will 
have the same structure as in the unmagnetized case, 

(8.80) 

One can then follow exactly the discussion in Sec. 14.21 to arrive at the integral 
equation (|4.^M|I for !Fl{x). Solving the latter looks a rather formidable task, given 
that Tk itself has the relatively complicated form H8.78|l . Remarkably, however, the 
solution can be found in closed form and is given simply by 

^..(x) = ^.(^-'^)/2 + |_^,-./2 (8.81) 

We were led to this result initially by a systematic series expansion of both Tk {x) and 
Tl{x) in terms of (x — 1)/ x. We do not detail this here, but verify in [Appendix B| by 
direct calculation that (|8.81|) does indeed solve (|4.34() . 

We next calculate the kernel M{t,t^). One inserts the scaling l|8.8U|l into (|8.()0|) . 
subtracting off and adding back on the contribution from the equilibrium part of L^^h 

m-\t)M{t,t^) = 5{t-t^) + 2T,- f dt' L^^^\t' - t^) 



+ / dt' L^J'^it' -t^) 



(8.82) 



This is done to account for the leading order cancellation of the second and third 
terms: 

/t poo \ 

dt'L^-'^it'-t^) = 2Te-L(2)(0)+ dt'Lf^{t'-t^) = ^(i_t„)('^-4)/2(8.83) 

where in the last step we have restricted ourselves to the aging regime t — <w ^ 1 
and used the asymptotic behaviour 14.23|l . ioq^(t' — t^) ~ \d{t' — tw)^'^~^^^^- In the 
remaining integral in H8.82|l . the factor 1 — J^L{t' /t^) ^ [t' — ty^)/tvr ensures that the 
integration no longer has its weight concentrated near t' — t^; put differently, after 
scaling the integration variable by to y = t'/t^j the integral is convergent at the 
lower end. We thus obtain in the aging regime 

m-\t)M{t, t^) = d{t - t^) + A,4''-4)/2 {-^(x - l)('^-4)/2 



+ / dy{y-lY'-'y'[l-TL{y)]\ (8.84) 



Inserting l|8.81|l for Tl, the integral in (|8.84|) can be done explicitly to give [2/(4 — 
^-jj2;(2-'i)/2(2; x)(<i-4)/2 ^Yic sum of the terms in curly brackets. A little care is 
needed here because the separate integrals over (y — l)('^~6)/2 l)(''-6)/2jr^j-y-j 
are divergent at the lower end. One can avoid this by analytical continuation from 
d > 4, where these divergences are absent, or by integrating from 1 + e to a; and 
taking e at the end. The 5-term is again subleading for long times in the relevant 
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combination (18.621). which we can write as 



5{t - 1^) - M(t, tMt^) - s{t - 1„) - / , ^f/;,, - ^Y'-'^^' 

= 6{t-t^)-^J^M{x) (8.85) 

with 

= (8.86) 

Here we have ehminated the constants Ad and fid using the foUowing argument. 
From l|0) . g{s)/m'^{0) = Zcq(s)/s^ for smaU s, i.e. Lcq(s) = s25(s)/m2(0). In 
the time domain this gives at long times L^{t) — —Lcq{t) ~ -~{dtf'm~'^{t) = 
-,i-\dt)H^^~^y^ = -Ai^i[(d-2)/2][(d-4)/2]t('^-6)/2^ go that AdMd - (d-2)(4-d)/4. 

Reassuringly, for d — > 4 the result (|8.86|l tends to Tm{x) — x^^/^, matching 
smoothly onto the result (|8.62() we found earlier in d > 4. For d < 4, the scaling 
function diverges as x — > 1. Since from l|8.60(l the continuous part of M{t,tv,) is 
exactly given by 2Tcm{t) for equal times, this indicates that the above aging regime 
expression must break down eventually when t — t^i becomes small, as expected. In 
the integrals where M appears below such effects can be neglected, however, because 
they only give subleading corrections. 

With the scaling of in hand we can now compute the connected correlation 
function Cm(t, iw) including non-Gaussian corrections, given by H8.29|l . After rescaling 
the integration variables the first part H8.31|l can be written as 

CW(i,tw)= / dy[5{y~x)-x-^TM{x/y)] 



dy^ [S{y^ - 1) - .FAf(l/j/w)]Co(iw2/,iwy') (8.87) 
where the Gaussian magnetization correlator is 

Coit^y, t^y') = {iT,t^/d) min {y{y/y^)^''~^^^\ yAy^/vY"'^''^^} (8-88) 
from H8.15|l . At first sight 18.87|l suggests, e.g. from the 5{y — x) term, an asymptotic 
decay of Cm"* ~ iwa^*^^"''-'^'' for large x. But this would not match continuously with 
the result iPTTHIl we found for d > 4. A cancellation of such leading order terms must 
therefore occur for a; — > oo. To show this explicitly, we verify from H8.86|l the identity 

■ dy [S{y ~x)- x-'Tm (x/y)] y-^^-'^^ = 

2a; Jo 

= ^(2-d)/4 _ ^^(2-d)/4 _ ^)(d-4)/2 ^ (8_9o) 

2 Jo 

Multiplying this by [ATct^j / d) dy^j [5{y^x ^ 1) ~ J^M{^/yw)\yv!'^'^^^'^ and subtracting 
from (|8.87|l exactly cancels all contributions in the range y > y^,, giving 

4T t Z"^" 
^w) = I dy„ I dy Tm {x/y) [5{y^ - 1) - Tm (1/yw)] 



°^ ,(.-2)/4 . x(.-2)/4- 



(8.91) 
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This shows that Cm'' /t^r is a scahng function of a: = t/t^. Its fuU a;-dependence 
has to be found numericaUy from (|8.91|) or via series expansion |22| . but we can 
obtain the large- cc behaviour that is required for the asymptotic FDR X°° in closed 
form. For x ^ oo one can replace the function J^M{x/y) with its asymptotic form 
[(d- 2)/2](x/y)(2-'')/4 from to get 



2(d- 2) 



-(d+2)/4 



T f 



, 2/w 
Vw — 





(d-2)/4 



y — 



(d-2)/A 



5{y^ - 1) - Tm 



2(d- 2) 
d + 2 



-(d+2)/4 



1 - 



2 r(d+i) 



(8.92) 



(8.93) 



This exhibits the expected leading order cancellation for large x, which gives an 
additional factor of 1/x compared to the naive result a;*^^"'*'/'*. 

To complete the calculation of the correlation function we need to evaluate Cm 
from (|8.32|) . which cannot be neglected for d < 4. This requires the long-time 
behaviour of CC{t,t^), which is given by (|7.5|l and (|7.()|) for t > t^ and t < tw, 
respectively, as for the unmagnetized case. The only modification arises from the 
different behaviour of g{t). One thus finds 



CCit,t^) 
CC{t,t) 



Git/t^), g{x) 



J dw U;('^-6)/2jr2 (^)g-2(x-l)u, 

a;/dww('^-6)/2jc-2(^) 
a;('^-6)/2g(l/a;) 



for a; > 1 
for a; < 1 



(8.94) 



The scaling of the equal-time value of CC is, from (|73|l . CC{t,t) = 7dt(''"'''/^ with 
Id = f^dT^ J dw w^''-~^^^'^Tl,{w); compare (|A.1|I . Inserting this into H8.32|) gives 

^ ^ dt' dt'^ M{t, t')M{t^,t'JCC{t', t'J 



1 



-idt 



' m(t')m{t' ) 7 -^^-f ~ ) C'C' (i' , t' ) ^ O 



^ f' dy^yi'-^y^TMilM r dyTM{x/y)y^^~^^'^g{y/y 

^l^dX Jo Jo 



c^yw2/w^A/(i/yw) 



du.FM(x/uyw)w(«-'^)/^-^^(u) 

2Md 



(8.95) 
(8.96) 
(8.97) 
(8.98) 



In the second line we have used 18.85|l to write M(t, t') = m{t)d{t — t') — 
t~^m~^{t')J^M{t/t') up to negligible corrections, and then immediately discarded the 
(5-function contributions, which give subleading corrections. 

Let us denote by U the value of the u-integral in H8.98|l . Since g{u) is defined 
separately for u > 1 and m < 1 in (|8.94|l . one splits the integral accordingly: 



2ndU 
Id 



x/y^ 



= / duJ^Mix/uy^)u^''-''^/''g{l/u)+ / duJ^M{x/uy^)u^''-'''>/*giu) (8.99) 



dw.FA/(a:V2/w)u-('^+2)/4g(u) 



x/y^ 



duTM{x/uy^)u^'^-'^^^'^g{u)i8.100) 
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We now need G{u) for u > 1. The denominator in (|8.94|) is Jd/ {(^dT^), and bearing in 
mind the definition H8.f 0|l of Tc with a = (c? — 2)/2 gives 



Id 



d/2 



(8.101) 



(8.102) 



(8.103) 



We now insert this into (|8. 100(1 and simphfy the numerical prefactors by using 



= ^(d — 2)(4 — d)/4 and the exphcit expression (|A.2|I for Xd to get 



U = 



Tr. 



r(V)r(M) 



duTM{xu/y^)u''^''+^^^^ 



X dy dy'{yy')(^~'y^l-y-y' + uy 



-d/2 



duTMix/uy^)u^^-''y^ C dy C dy' 



(8.104) 



This is the u-integral from H8.98|l and so overaU we have a four-dimensional integral 

f 2) 

over y^,u,y,y' for Cm . In general this cannot be evaluated in closed form; a series 
expansion is given in |22| . The large-x behaviour, which will give us the asymptotic 
FDR, is easier to extract. In the first it- integral of H8.104|) one can directly use the 
asymptotic form of !FM{xu/y„). One can show that for large x the same replacement 
can be made in the second integral, and the upper integration limit sent to infinity 
thereafter. This gives for the large- a; behaviour of U 
d - 2 T,Vd 



U = 



2 r(^)r(^) 



(2-d)/4 (d-2)/4 
i/w 



where Vd is a d-dependent numerical constant given by 



Vd 



du{u-^^+^^'^ + 1) dy dy'{yy'f''-^'^'\l^y-y' + ur'"^ 



(8.105) 



(8.106) 



Inserting (|8.105|) into 



d-2 



the remaining j/w-integral can be done explicitly to give 

T{'^)Vd 



T{d+l)T{^) 



Tt T-^'^+'^^Z^ 



As anticipated this has the same scaling as the first contribution 
correlation function, so that overall for large a; 



(8.107) 
to the 



C (t t )- C(i) + C(2) - ^ t a-"('^+2)/4 



d + 2 



1 - 



d-2r{^)T{^) 



r(d + i) 



d-2 r{^)Vd 
2 r{d + i)r{^) 



(8.108) 



The magnetization response function is rather easier to find, by using (|8.14|l 
and H8.85|l in H8.39|l and rescaling the integration variable to y — t' /t„ as usual: 



Emit, tw) = / dy [Siy - x) ~ x-^TMix/y)] y-^'^-^^^ 



(8.109) 
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The structure of this is rather similar to Cm ^ , and by subtracting the vanishing 
term (|8.9U|) one again gets a significant canceUation, 



Inserting the exphcit form 



of J^]\i one then gets simply 



1 



(d-2}/2 



which for x oo behaves as 



d-2 



(8.110) 



.111) 



.112) 



With the expressions (|8. 108(1 and ((8.112|l for correlation and response in the 
limit of long, well-separated times we can now finally compute the asymptotic FDR 
defined by 



X:, 



lim 



^c-^m(^; tw , 

C^{t, iw) 



d_2 r{^)v. 



d + 6 



-I -1 



d + 2 



1 



2r(^)r(^) 
r(d + i) 



i.ll3) 



2 r(d + i)r(^) 

where Vd is given by (|8.10t)|l and the prefactor 4/((i + 6) accounts for the fact that 

an ~ iw2:~(2+d)/4 _ t(^rf+6)/4 ^^^^^ C^{t,t^) = [{d + 6) / 4t^]a^{t , t^) . 

Before exploring this result, let us comment briefiy on the limit as rf ^ 4, which 
should make contact with our results in the previous subsection. The contribution 
from Cm , which appears in the second line of (|8.108|l and (|8.113|l , vanishes linearly in 
4— d in this limit because of the factor r~^{{A—d)/2); Vd stays finite as we show below. 
This is consistent with the fact that in dimension d > 4 this term does not contribute. 
In the limit c? — > 4 one has, from H8.108|l . Cm = TctvtX~^^'^ /2 for large x which matches 
precisely (|8.70|l for d > A. Similarly, the large- a; magnetization response (I8.111|l for 
d ^ 4 is i?m = x^^^'^ in agreement with H8.72|l . 

We now look in more detail at the d-dependence of the asymptotic FDR (|8.113|l 
for the magnetization. Expanding in e = 4 — d one has 

-1 



4 / 28 ^4 

5 ^ V225 ~ y 
where V4 is the limiting value of Vd for d 



+ 4- 
/ 28 
V225 



2 



+ -Vi + 0{e') 



.114) 



Vi = -3ln2 



11 
12 



TT 

24 



so that 

4 

^m -5 



21n2 
15 



53 
900 



nnn 120 ' ^ ^ 



which can be worked out explicitly as 

(8.115) 

(8.116) 



It is remarkable that in a system as simple as the spherical model, where standard 
critical exponents are rational functions of the dimension d, the magnetization FDR 
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Figure 4. Asymptotic FDR for the magnetization, for critical coarsening 
with nonzero initial magnetization. Solid line: Full theory I8.113i including non- 
Gaussian corrections; dotted lines indicate the first-order expansions near d = 2 
and 4. Dashed line: Gaussian theory J^H}. 



for magnetized initial states is very much more complicated, and irrational already to 
first order in e. 

In the opposite limit d ^ 2, the u-integral in the definition (|8.1U6|I of Vd diverges 
at the upper end; dropping all non-divergent corrections gives the leading divergence 
as 

du dyj^ dy'il-y-y' + u)-"/^^ duu-"'^ = j-^ (8.117) 

This divergence balances the vanishing prefactor (d — 2)/2 in (|8.108() and (|8.113|l 
whereas the contribution in the round brackets coming from Cm'' vanishes linearly 
with d — 2. Consequently the asymptotic behaviour of the correlation function is, for 
d close to 2, Cm = [{d ~ 2)/2]Tciwa;-^ and the FDR becomes = 4/(d 6) = 1/2. 
One can also obtain the leading order correction in e' = rf — 2, which is given by 

^m + + (8.118) 

The only subtlety here is working out the subleading term of Vd- Setting Vd — 
2/e' + flo + . . ., flo can be obtained as the limit for c? ^ 2 of H8.1()6|l with u^'^/'^ 
subtracted from the integrand. The limit can be taken in the integrand itself for 
finite u, giving a convergent integral with value 3/2 — 7r^/12. But one has to account 
separately for the large-u tail 1 + / dy dy' (yy')^'^^^^^^] which integrates to 

{-l + A/d^)[2/{d-2)] -2 for d ^ 2, giving oq = -1/2- 7rVl2 overall. 

In summary, the asymptotic magnetization FDR X^ decays from 4/5 in d = 4 to 
1/2 in ci = 2. Fig. 01 also shows numerical values for intermediate d. X^ is larger than 
the FDR (|8.12|) from the Gaussian theory except in the limit d ^ 2; it also remains 
larger throughout than the FDR ()3.3|) for unmagnetized initial states, shown in Fig.|21 
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Figure 5. Normalized magnetization FD plot for d = 2,3, 4, showing normalized 
susceptibility Xm versus normalized correlation Cm in the limit of long times. For 
d = 4 the plot is a straight line with (negative) slope 4/5 as expected. Increasing 
deviations from this appear as d decreases towards 2. 



We next turn to the shape of the FD plot for the magnetization. As explained in 
the introduction, this is obtained by plotting the normalized susceptibility Xmit, iw) = 
TcXm{t,t^)/Cn,it,t) versus C^{t,t^) = C^{t,t^) / C^{t,t). In the limit d 4 the FD 
plot must be a straight line with slope = 4/5, by continuity with the results for 
d > 4. Numerical evaluation (see Fig. shows that, as the dimensionality decreases, 
the FD plots deviate progressively from this straight lines. The most extreme case is 
the limit d — > 2, where analytical forms can be found. 

To find the correlation function for c? 2, it is useful to note that the 
scaling function !Fm{x) becomes equal to 5{x — 1) in the limit. Formally, one sees 
easily from (|8.86|l that for any smooth bounded function f{x) and fixed c > 1, 
dx J-M{x)f{x) /(I) because the divergence of J^]\i{x) at a; = 1 becomes non- 
integrable in d = 2. We now exploit this to simplify U from H8.104|l . In the first term 
in square brackets, the argument xu/y^j of J-^m is always larger than x/j/w and hence 
than 1 (except at the irrelevant boundary = x of the j/w-integral); in the limit 
d — > 2, this contribution therefore vanishes. In the second term, replacing J-Mix/uy^) 
by S{x/uyv, — 1) and expanding the prefactor to leading order in (d — 2)/2 gives 



d~2 



U_ 

Performing the integrals over y and y' , one has 



dy / dy' {1-y-y' + x/y„y 



U_ 

T'c 



d~2 



In 1 - 



In 



(8.119) 



.120) 



(2) 

We can now use this limit form of U to get the contribution Cm to the correlation 
function; recall that U was defined as the u-integral in H8.98|l . In the remaining 
integral in this equation one can again replace J^m[^IVw) by 5{l/y^ — 1) so that only 
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1 contributes, giving 



d-2 



a; In 1 



In 



.121) 



One can show that the other contribution Cm'' to the magnetization correlator, given 
by (|8.91|l . vanishes as ~ (d — 2)^ for d ^ 2 so that to leading order Cm = Cm'' 
as anticipated in writing H8.121() . The corresponding response function is found by 
expanding H8. 111(1 to leading order in d — 2: 



Rn 



d-2 



In 



1 



To get the FDR it only remains to work out the t„-derivative of C„ 



Cm = dtCm — 



d-2 



■TAu 



1 



1 



We therefore get, for the full dependence of the limiting FDR for d 
X — 1 1 



= In 



1 



In 



1 



1 



(8.122) 

(8.123) 
2 on scaled time 

(8.124) 



For a; ^ oo this gives — \/2 consistent with the discussion above. In the limit 
a; — !■ 1 of comparable times, on the other hand, Xa\{x) approaches 1, logarithmically 
slowly; the FD plot for d ^ 2 therefore starts off with a pseudo-equilibrium slope. 
Interestingly, this implies that the trends of the slope with d are different at the two 
ends of the plot: for well-separated times {x oo), the slope decreases from 4/5 to 
1/2 as d decreases from 4 to 2; for comparable times {x —^ 1) it increases from 4/5 to 
1. 

To get the FD plot itself we need the susceptibility, which is found by integration 
of the response (18.1111) as 

ft 



Xm(i,iw) 



dt' R„,{t/t') =t dzRn,{l/z) 



1/x 



1/x 

(1 



.(d-2)/2 



Expanding to linear order in d — 2 and integrating gives 



Xmit,t^ 



2 X 
— t — 



1 



1 - In 



1 



.125) 



(8.126) 



.127) 



The normalized susceptibility Xm is obtained by dividing by Cm(i, which 
from H8.121|) for x 1 equals 21n2[(d - 2)/2]i. The resulting FD plot for d ^ 2 
is shown in Fig. \S\ together with the ones for d — 3 (determined numerically) and 
d = A. For d^2 the approach of the slope to the equilibrium value X = 1 for x — > 1 
is difficult to see because it is logarithmically slow. As expected from the trends with 
d in the initial and final slopes of the FD plot, the curve for d ^ 2 is the most strongly 
curved, while in d = 4 we have the anticipated straight line with slope 4/5 that is 
required by continuity with the results for d > 4. 

To quantify the shape of the FD plot further one can also consider the axis ratio 
Y, defined as the limiting value of Xm(i, ^w) = 7cXm(i, iw)/Cm(i, t) for well-separated 
times i » iw S> 1. In equilibrium this would correspond to the FDT for the static 
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quantities, i.e. equal-time fluctuations and static susceptibilities. Out of equilibrium, 
if the FD plot is straight then Y coincides with X; if it is not, then is has been 
argued that in some circumstances Y can be more relevant for characterizing effective 
temperatures than X |2Sj- For the magnetization FD plot in the current scenario of 
magnetized initial states, we see from Fig. that Y decreases along with X^ from 
4/5 as the dimension is lowered below d = A. The two quantities begin to differ more 
noticeably as d decreases further, with X^ = 1/2 and Y — l/(21n2) = 0.7213. . . in 
the limit d ^ 2. 

9. Summary and discussion 

In this paper we have considered the non-equilibrium dynamics of the spherical 
ferromagnet after a quench to the critical temperature T^. Our focus has been 
the calculation of correlation and response functions and the associated fluctuation- 
dissipation ratios (FDRs) X(t,tvi)- The key quantity that can be extracted from the 
latter is the asymptotic FDR X°° for large and well-separated times t ^ 3> 1; it is 
independent of model-specific details within a given dynamical universality class. We 
were motivated by two questions: how does X°^' depend on the observable considered, 
both with regard to the lengthscale and the type of observable (spin, bond, spin 
product)? And what is the effect of initial conditions, in particular the presence 
of a nonzero magnetization in the initial state? The first question has implications 
for the interpretation of T/X°° as an effective temperature, which is plausible only 
if this quantity is observable- independent. The second one allowed us to uncover 
whether different initial conditions can lead to different universality classes of critical 
coarsening. 

A peculiarity of the spherical model is the weak infinite-range interaction 
produced by the spherical constraint. This requires that one distinguishes between 
long-range or "block" observables, which probe lengthscales large compared to the 
(time-dependent) correlation length but small compared to the system size, and 
global observables whose behaviour depends on correlations across the entire system. 
Technically, the first case is much easier to treat because the standard theory where the 
spins have Gaussian statistics can be used. Global correlation and response functions, 
on the other hand, require non-Gaussian corrections arising from the fluctuations of 
the effective Lagrange multiplier. 

We dealt with the case of flnite-range (i.e. either local or long-range) observables 
in Sec. 13 For spin observables, we found in the long-range case H3.3|l the same 
X°° as for local spin correlations and response [H]. This was as expected from the 
general correspondence between local and long-range observables discussed in the 
introduction. The FD plot for the long-range spin observable, i.e. the magnetization, 
is a straight line in the long-time limit. This is as in the Ising case in d = 1 ^0], but it 
is interesting to note that here it holds for all dimensions d > 2. On the other hand, in 
the Ising model with d>2, KG arguments have been adduced ^lElEl to suggest 
that the magnetization FD plot should not be straight, though with deviations that 
are likely too small to be detectable numerically ^Sl- It is likely that the Gaussian 
statistics of the spherical model are responsible for producing a simpler, straight-line 
magnetization FD plot in all dimensions, although it would be interesting to know 
whether any other models have this property. 

We then looked at the effect of the type of observable on X°°, considering both 
bond and spin product observables in either the local or long-range versions. The 
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results in equations H3.15|l . (|3.2Ufl . ()3.2tifl . H3.31|l show that, although the precise time- 
dependence of X{t,t^) varies, the asymptotic FDR X°° is the same in all cases. This 
is consistent with general arguments J3| suggesting that for a Gaussian theory all 
observables should yield the same X°°. In contrast to the Ising case jTO], not all long- 
range observables give nontrivial FD plots; in fact, only the block product observable 
does so, and only for d < 4, while all others produce pseudo-equilibrium FD plots for 
long times. 

The bulk of the paper was concerned with the more challenging analysis of global 
correlation and response functions, focussing mostly on the energy as a key observable. 
In Sec.01wc constructed a framework for calculating non-Gaussian corrections to the 
spins, which are 0{N~^^^) to leading order. This lead to the general expression (|4.11ll 
for these leading-order corrections. It involves a two-time kernel L{t,t^) which 
from H4.9|l is the functional inverse of K{t,tvr) defined in (|4.8|l . The basis of all 
subsequent calculations is the determination of the long-time scaling of these two 
functions, as summarized at the end of Sec. 14.21 

In Sec. 13 we obtained general expressions for energy correlation and response 
functions, in terms of the kernel L and other quantities known from the Gaussian 
theory; the results can be found in eqs. H5. 515. 1215. 1415. 19ll and (j5.29(l . Evaluating 
these first for the equilibrium case, we found that the energy correlation and 
susceptibility display a plateau for T just below Tc and d > A; this is caused by 
the q = wavevector, i.e. by the slow relaxation of the global magnetization. In 
Sec. El we proceeded to the long-time analysis of energy FD behaviour in the non- 
equilibrium case for d > 4; the key results are (|6.23|l and (|6.24(l . The associated 
FDR is given explicitly in l|6.22|l and has the same asymptotic value X°° — 1/2 as 
for all other (finite-range) observables in ci > 4. The analysis of the case d < 4 is 
more difficult, and we were able to find closed- form results H7.15I7.16|I only in the 
limit of well-separated times t/tv, ^ 1. This is, however, enough to determine X°°, 
with the result H7.19(l . Evaluating this, both numerically and by expansion in 4 — d 
and d — 2, the crucial conclusion is that it does not coincide with the asymptotic 
FDR for finite-range observables; see Fig. O A naive interpretation of T/X°° as 
an effective temperature for critical coarsening dynamics is therefore ruled out, since 
such a temperaure ought to be observable-independent. On the other hand, to first 
order in 4 — d the result agrees with an RG calculation jEl for the 0(n)-model. We 
conclude that non-Gaussian corrections to the FD behaviour of global observables in 
the spherical model capture genuine physical effects that have close counterparts in 
more realistic systems with only short-range interactions. 

Finally, in Sec. |H1 we turned our attention to critical coarsening starting from 
magnetized initial states; physically this situation could be produced by an up-quench 
from an equilibrated state at a starting temperature T < T^. We concentrated on the 
simpler spin observables and found that already for them, the presence of a nonzero 
magnetization makes global properties sensitive to non-Gaussian corrections. As with 
the energy fiuctuations, it is the global correlation and response functions that make 
contact with the results for short-range models, as obtained recently for the Ising 
case 1201: we find X"^ = 4/5 for d > 4, eq. (|8.73(l . This is distinct from the value 

= 1/2 for the unmagnetized case, indicating that magnetized critical coarsening 
is in a separate dynamical universality class. Surprisingly, even the expressions for 
correlation and response functions themselves, which are not expected to be universal, 
coincide with those for the Ising case. It remains to be understood whether this is 
accidental or has more profound origins. For the case d < 4, we obtained new exact 
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values for the asymptotic FDR of magnetized critical coarsening. The magnetization 
response (18.11111 can be found explicitly for long times, while for the magnetization 
correlator only the asymptotics for well-separated times H8.108|l can be written in 
closed form. The resulting X^, eq. (|8.113|l . is surprisingly nontrivial: while it matches 
continuously with = 4/5ind>4 and approaches the simple value X^ — 1/2 
for d — > 2 as shown in Fig. 01 it is irrational already to first order in an expansion in 
4 - d, eq. dHUa, or d - 2, eq. ^im. 

While the conclusion of our calculation as regards the existence of a well- 
defined effective temperature for critical coarsening is negative, the issue of dynamic 
universality classes and new asymptotic FDRs due to magnetized (and possibly other, 
different) initial conditions clearly deserves further study. Results for systems with 
short-range interactions, such as the 0(n) and n-vector models, would be particularly 
welcome. After the present work was completed we became aware that a first step 
in this direction has recently been taken by the authors of Ref. [23 , who calculated 
the FDR for the n-vector model with a magnetized initial state within an e-expansion 
around d — 2. Intriguingly, their result X°° = 1/2 for d = 2 itself agrees with ours, 
but the first-order correction in o? — 2 remains rational even for n — > oo. It therefore 
disagrees with our spherical model result H8.118|l . This appears to be the first example 
of genuine differences between the spherical and n-vector (with n — > oo) models, which 
are known to have identical properties in equilibrium |25| and within a Gaussian theory 
of the dynamics. 

As regards future work, we note first that a complete classification of dynamical 
universality classes within critical coarsening remains to be achieved. An earlier study 
of the spherical model considered initial conditions with long-range correlations but 
no overall magnetization; this yields no new (non-zero) values of the asymptotic FDR 

|26| . The presence of a non-zero magnetization thus appears to be important 
for observing new phenomena, and is reflected in our calculation by the fact that 
non-Gaussian fluctuations become important. Whether there are yet other initial 
conditions that could give rise to distinct values of X°° is an open problem. 

Our general framework for treating non-Gaussian corrections to the dynamics 
can also be applied in other contexts. For example, it can be used to analyse 
the fluctuations across thermal histories of correlation and response functions 
that have been coarse-grained across a finite-sized system. The properties of 
these fluctuations should be useful for understanding dynamical heterogeneities in 
coarsening dynamics j27| . and we will report on the results of such a study shortly. 
We have also extended our approach to non-Gaussian corrections for the dynamics 
of e.g. the 0(n)-model with large but finite n, opening up the attractive prospect of 
obtaining exact results analogous to the ones in this paper for models with exclusively 
short-range interactions. 
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acknowledges the hospitality of the KITP, Santa Barbara, where this research was 
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Appendix A. Evaluation of X^ 

In this appendix we evaluate the various numerical factors in the asymptotic 
FDR H7.18|l for the energy in d < 4. For (3d we already have an expression (|7.14() . 
By definition, 7^ is the large-i hmit of CC{t,t)/t^^^'^^^'^ . Its value can be deduced 
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from H7.3|l : for small q, q — lu 
proportionality constant, 
because it cancels from the overall result 



^ and hence (dq) — a^dujuj'''^ (The value of the 

(47r)^''/^r^^((i/2), is not actually needed expHcitly 
Thus from ((73|l 



Next, Ad is defined by Liq^(i) = Xdt^'^-^'^/^ for large t. 



?(2)/ 



Id = <TdT^ jdww'^^-'^'^TUw) (A.l) 

The Laplace transform ot 

L^cq therefore has a singular term Lcq^s) = Xd^{{d — 4)/2)s^^~'^)/^ for s 0. 
From H4.21() . this must match the corresponding singularity in —1/Kaii{s). The 
latter follows from KTI^ as i^cq(s) = s^'^-'^)/^ j.^^^ j^^ ^(d-4)/2(;^ _^ 2w)-i = 
g(d-4)/22.^^^2(2-d)/2p((4„ ^)/2)r((d- 2)/2). This tells us that 

= -r,a,2(2-'^)/2r (^) r (^) r (^) (A.2) 

Finally, we need ad from H7.13|l . Inserting (|2.25|l into H7.7() and using (|A.ip . the 
scaling function Q{x) can be written for x > 1 as 



g{x) 



Id 

r(|)2(4-rf)/v^r,2 
Id 



dyy 



dy 



(d-4)/2 



dy'{yy'Y''-^^^^e-^'^'+''-y-y'^^ 



y/2 



d/2 



_ 2 r(f)2(4-'^)/2^^r,^ 

^ - 2 7rf 

We define the integral as G{x) 
5(1) = 1 one then has G{x 
ad, ad = adG{l), is 

-d/2 



(d-4)/2 



_ ,A(2-rf)/2 



(x-y) 



(A.3) 
(A.4) 

(A.5) 



dyy , , 

1 + a; — y 

It is an unnormalized version of G{x)] because of 
G{x)/G{i)- The corresponding unnormalized value of 



ad 



dx X 



^9d 



9dX 



dx xG (x) 



G{x) 



dx 



-gdx-"'^ - G{x) 



(4-d)(rf-2) 



dxx'^-'+^y^Gix) 



dx X 



Ix-'^l'^ 
d~~2 



G{x) 



(A.6) 
(A.7) 
(A.8) 



Here we have used (I7.7|) to express the values of G{x) for x < 1 in terms of those for 
a:; > 1; we also inserted gd — gdGi^) — i^nix-,ooGix)x'''^^ — 2/id— 2) which follows 
from (|A.5|I . The first integral has no divergences. The second one could in principle be 
left as it is for numerical evaluation, but it is useful to rewrite it using a dimensional 
regularization trick, as follows. If we extend the definition of ^(a:) to d > 4 using ljA.5|) . 
then both parts of the second integral in (|A.8|I are separately finite for 4 < d < 6, 
and we can evaluate them first in that range of d and then analytically continue to 
d < 4. (The region a; « 1 causes no difficulty since ^(a;) ~ (a; — l)(4-'i)/2 fQj. ^ ^ 1, 
which remains integrable for d < 6.) The first part gives [2/{d — 2)][2/{d — 4)] and 
just cancels the first term in l)A.8(l . so that 



dx x 



2x-'^/2 



d 
dx 



Gix) 



{4:-d){d-2) 



dyy^^-^y\x-y)^^-^y^- 



sin[7r(d-4)/2] 



+ / dx dyy 



,{d-4)/2 



x-y 

(a;_y)(2-rf)/2(l 



y) 



1 



y 



(A.9) 

(A.IO) 
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The cancellation of the pure power-law term proportional to g^j from the integral (|7.13|l 
is a feature well-known from dimensional regularization in field theory. Inserting the 
last expression into IjA.Sp then gives H7.20|) in the main text. 

Collecting the above results, the asymptotic FDR H7.18|l for the energy in d < 4 

is 



dad 



2 r (f ) 2(4-<i)/2a,T2 / r {^) r {^) 



d-2 7rf \^ r(d) 



2 r(d)r(i^) 

dad r {^) 
which is the result stated in H7.19|l . 



Id 



(A.ll) 
(A.12) 



Appendix B. Solution for J-l{x) for magnetized case in d < 4 

In this appendix we prove that the solution J^l{x) = 2/(4 - d)x'-^^'^''^^ + (2 - d)/{4: - 
d)x^'^^^ given in l|8.81|l does indeed satisfy the integral equation H4.34|l 



dy {x - vY^~^^'\v ~ li'^-^y^ [Tk [x/y) {y) - Tk{x)] = 



(B.l) 



with d < A and Tk{x) given by (|8.78|l for the magnetized case considered here. 

The difference in the square brackets can be written as the integral of a partial 
derivative 

(•V 



(x/y) Tl [y) - Tk{x) 



dzdz [Tk {xIz)Tl [z)] 



(B.2) 



so that, after exchanging the order of the integrals, the left-hand side of l|B.ip becomes 



I = 



dzdA^K{x/z)TL{z)] / d2/(x-y)(2-'')/2(y-i)('i-6)/2 



1 



' dz{x^ z)(4-'^)/2(z - l)('^-4)/2a, [Tk{x/z)Tl {z)] 



(B.3) 



(B.4) 



x- 14-d, 

where in the second step the y-integral has been performed. Eq. HB.4|I can equivalently 
be obtained from (jB.l|) by integrating by parts. From (|8.78|l for Tk{x) and the form 
given above for Tl{x), we can find the z-derivative of Tk {x/ z) Tl{z) explicitly. After 
simplifying the latter, equation ljB.4p becomes 



l = r, r dz{x- zY^-^y\z - l)('*-4)/2 /_i + ^ + 
./1 \ z 2z^ 



d-2 



—X — z{x + 1^ 



X Z-('^+2)/2(a; - z)(^-4)/2 / _ y^(2-d}/2yid-4)/2 

J z/x 



where 



^(2-d)/2 2(ri„2) 
x-1 (4-(i)2 



(B.5) 



(B.6) 
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collects all the prefactors. By rescaling y by a factor x, the last integral in HB.5() can 
be transformed to dy {x — i/)(2-<i)/2y(d-4)/2 interchanging the y and z-integrals 
then gives 
I 



dz {x - zY^-'^y^iz - (^-^ + ^) 

dy(a;_y)(2-<i)/2j^(<i-4)/2 



X I' dz - ^ - Zix + 1) 

The integral on the last line can now be calculated explicitly to give 



giy) 



1 



y 



(d-2)/2 



d-2x 



2 y 



1 



Relabelling y to z, eq. IjB.Tp can thus be written as a single integral 
I 



dz (x-zY^-''^^ 
dz (x-z)(2-'i)/2(^_i)(rf-4)/2 



{z-ir-'y\x-z) f-l + AA+z^'^-^y^giz) 



X 
72 



d-2 A-d 



(B.7) 
(B.8) 

(B.9) 
(B.IO) 



Finally, the variable change v = {x/ z — l)/(a; — 1) transforms the integration range to 
. . . 1 and leads after a few simplifications to straightforward Beta-function integrals: 

^^'^ (B.ll) 

(B.12) 



I = x(2-'')/2(^ _ 1) dv V^'-^'^' + „^ (1 _ ^,)('i-4)/2 



- rr('2~d)/2 



(x-l) 



4-rfr(^)r(^_)_^r(^)r(^) 



2 r(i) r(2) 

This proves that IjB.ip is indeed satisfied, as required. 
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